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Section I 

Introduction 


The development of mathematics in the nineteenth century was 
so extensive that few mathematicians after Gauss could hope to 
understand more than a small sector of the subject. Outstanding 
in the mid nineteenth century was the development of different 
algebras. For centuries algebra had been applied to numbers: 
Descartes and Newton had developed the algebra of points 
represented by numbers; but matrix algebra is the algebra of 
coefficients of linear equations, such linear equations representing 
operations or transformations. 

Most of the symbols used in matrix algebra are the same as 
those used in number algebra, but the operations they represent are 
different or the quantities represented are not the simple numbers 
of common usage. Before proceeding to the explanation of what 
matrices are and before defining the symbols used, it is necessary 
to be clear about the fundamental laws of number algebra. 

In the primary school the child learns from handling number 
apparatus that 3-1-4 is the same as 4 + 3 and that the order of 
addition is not important : the result is the same whichever way 
he does the problem. The numbers 3 and 4 can change places, i.e. 
commute (think of the American use of the word commuter 
for a worker who changes places night and morning between his 
place of work and place of living). We say, in addition of numbers, 
that the addition is commutative. Similarly the child soon learns 
that 3x4 is the same as 4x3, i.e. multiplication of numbers is 
commutative also. When handling three numbers such as 2, 3 and 5, 
the addition of 5 to the sum of 2 and 3 gives the same result as the 
addition of 2 to the sum of 3 and 5, that is (2+ 3)+ 5 = 2+ (3 + 5). 
Note that the brackets used in this last expression denote the 
order of the operations. If we examine the statement (2+ 3)+ 5 
= 2 + (3 + 5) we see that the 3 can be associated with the 2 
first or with the 5, but the sum is obviously the same. In the same 
way (2 x 3) = 6 and 6 x 5 = 30 or (3 x 5) = 15 and 15x2 = 30, 
therefore the 3 can be associated with the 2 first or the 5 first 
and the continued product of 2.3.5 is still 30: the order of multi¬ 
plication of numbers is not important. We say that addition and 
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multiplication are associative when we are dealing with the 
algebra of numbers, i.e. (2 x 3) x 5 = 2 x (3 x 5). 

I hese relations can be expressed using symbols to represent 
numbers: 


/ Commutative Laws 
(i) For addition 
a + h = b + a 

II Associative Laws 

(i) For addition 

(</ + /?) + <• = </ + (/> + c) 


(ii) For multiplication 
a.h = b.a. 


(ii) For multiplication 
(a.b)c = a(b'C) 


By using simple numbers we can illustrate the third important law 
called the Distributive Law. 

Ill Distributive Laws 

(i) a(b + c) = ab + ac 

(ii) (b+c)a = ba+ca 


In number algebra the operations involving unit number 1 and 
zero 0 are important. Early in his number work the pupil learns 
that a number when multiplied by 1 remains unaltered in value 
and if zero is added to any number it again remains unaltered. 
For any number represented by a. 

*/ + 0 = a 
and 

a . 1 — a 


Later he learns that if —2 is added to -f 2 the sum is zero, or if 
— 3 is added to +3 the sum is zero, generally */ + ( — «) = 0. 
Then when he encounters the simple equation x + a = 0 he 
knows that x = —a because ( — a) + a = 0, but since addition is 
commutative this is generally written 

</ + ( — u) = 0 . 

In multiplication a further difficulty is soon met, when he tries 
to solve the equation x.a = 1 he finds that a cannot be equal 
to 0. Provided a is not equal to zero, then x is called the 
reciprocal of a and is denoted by a~\ therefore a.a~ l = 1. 
The process of subtraction in number algebra can be approached 
in two ways. Thus, directly, we can say k 7 minus 3 is equal to 4*. 
i.e. 7 — 3 = 4, but this could be turned into a process of addition 
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by asking ‘what number added to 3 gives 7\ i.e. 7 = 3 + 4. Gener¬ 
ally if a — b = c we say that c is the remainder but then a = b + c 
and c is the number which must be added to b to give a . The 
process of division can also be dealt with indirectly as soon as the 
meaning of the reciprocal is appreciated. Thus 3 + 2 can be written 
4 = 3x^- = 3x2 l , i.e. division has become a process of multi¬ 
plication. Generally for any number a divided by b the result can 
be obtained by multiplying the number a by the reciprocal of h, i.e. 


The algebra of numbers can be reduced to two fundamental 
operations of addition ( + ) and multiplication ( x) or (.). 

These ideas will be used in matrix algebra and where the symbols 
of number algebra are appropriate, they will be used in matrix 
algebra, but often with extended meanings. 

Matrix algebra can be approached in different ways, for the pure 
mathematician it can be started as an abstract algebra of groups 
or sets of coefficients, but for future scientists, engineers or tech¬ 
nologists a more practical approach is often more rewarding and 
meaningful: such approaches can be started in studying simple 
electrical theory or elementary equations or linear equations, but 
all these different approaches lead ultimately in the same direction. 
With younger pupils methods of coding messages using matrices 
arc often an exciting way of introducing matrices. 
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(2 x 2) Matrices 


We start with a study of (2x2) matrices in some depth. Many 
courses on matrix algebra plunge the student into great difficulties 
before he is familiar with the terms, concepts and processes: fre¬ 
quently the study of (2 x 2) matrices is brief, or neglected, or ignored, 
whereas it can provide the student with a rewarding introduction 
to some of the main ideas of matrix algebra and indeed the applica¬ 
tions and study of more complicated matrices are only extensions 
of the applications and ideas which can be acquired so easily by 
dealing with elementary (2 x 2) forms. In a simple way the student 
can gain insight into the structure of a new algebra of wide 
application and receive an introduction to many new concepts such 
as vectors. Any modern course of mathematics must deal with 
set theory and matrix algebra as its two main branches. 

The main practical applications of matrix algebra to linear 
programming, quantum theory, the theory of relativity, nuclear 
physics and above all to electrical network theory readily suggest 
some such practical approach to the subject. 

Meaning of a Matrix 

Consider the following arrangement. 


*1 


y i 

*2 

M 

yi 





Diagram 1 


Suppose M represents an electrical network with x, and x 2 
representing the inputs to the circuit and y, and y 2 the outputs 
from that circuit which is so constructed that y, and y 2 are each 
dependent on both x, and x 2 in a linear fashion. Such circuits can 
be constructed: this is the sort of condition that arises frequently 
in electronics, v, and x 2 can represent the e.m.f. applied to M or 
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the currents flowing into Af, similarly y, and y 2 then represent 
the e.m.f. or currents from the circuit in M. If V| depends on both 
x, and x 2 in linear fashion this can be expressed mathematically as 

y, = «,x,+/),x 2 where a,, b, are coefficients, i.e. numbers; 


similarly 

y 2 = o 2 x, +b 2 .\ 2 where u 2 . h 2 are coefficients, i.e. numbers 

These are linear equations of normal number algebra, on the left- 
hand side w'e have the outputs expressed as some operation on the 
inputs. Using a new notation we can express the relations differently 
in such a way that the two linear equations can be expressed in 
one matrix equation. If the two outputs y, and y 2 are put together 
yj I where the brackets have a special meaning in matrix algebra, 

then the coefficients a,, h, and a 2 , b 2 can be shown similarly as 


a, b; 

a 2 b 2 
matrix 


and the inputs x x and x 2 are put together in these special 


brackets 


*2 


Now the operation which transforms the inputs 


outputs 

>V 

is represented by the matrix 

a, V 


>’ 2 . 


a 2 b 2 


into the 


, the numbers 


a j, V tf 2 , b 2 are called the elements of the matrix. 

are called column vectors or column 


*1 

and 

y r 

.*2. 


>’ 2 . 


MATRICES. 

The. reason for the alternative use of column vector will be given 
later. 

The idea of representing the group of numbers in the correct 
order, which are responsible for transforming the inputs fX|l into 


the outputs 


in 


|*2 

a matrix form was due to cayley in the 


middle of last century. The numbers in the matrix are the co¬ 
efficients of .Yj and xj in their correct order they are ordered. 
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detached coefficients which have an algebra of their own 
called matrix algebra. 

Since this matrix a t ft,! has two rows and two columns it is 
a 2 b 2 

called a (2x2) matrix. In words—*a two by two* matrix. If this 
particular matrix is denoted by capital A and the two column 
matrices fx, 1, Vi ' by X and then the operation of the electrical 



1 

y i' 

.*2. 


>’2. 


Y= A.X 


where A represents an ordered set of numbers not one number 
and the product of matrix A with column matrix X must not be 
mistaken for ordinary multiplication, but a new kind of matrix 
multiplication (to be investigated later) which yields 

a i .x 1 +b t .x 2 

and 


a 2 .x l +b 2 .x 2 

Before dealing further with matrix multiplication, the process 
of addition and subtraction of (2 x 2) matrices will be explained 
or defined. Many processes are defined, not proved, the definitions 
giving required results. 

In future A will represent a matrix and the small letter a will 
represent a number in number algebra. 

Addition and Subtraction of Matrices 

Returning to the original electrical circuit: 




y 


M 

y 





Diagram 2 


Suppose that the circuit M is adjusted so that with the original 
inputs x, and x 2 , the outputs are now y' and y" such that 
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and 


y‘ = a,x,+fr,x 2 

y" = a 2 - v i +h 2 x 2 


In matrix notation 

y " 

= 

a \ hi" 

X,' 


y" 


cl 2 h 2 

X 2. 


Now suppose that the circuit M is altered so that an additional 
output y' x and yjj is secured in the output channels such that 

/. = a 3 .x,+b 3 .x 2 


>•; = a 4 -Xi+b*.x 2 


In matrix notation 

y'i 

= 

« 3 b 3 - 

x, 


ft. 


a 4 b A 

x 2 


or Y" = BX. 


If the new total outputs arc denoted by y, and y 2 as in diagram 1. 
then 

y, = y'+y, and y 2 = y"+y'i, 

y, = y'+y' x = a l .x i +b l .x 2 +a 3 x l +b 3 .x 2 \ 

= («i +a 3 )Xj +(b, + b 3 )x 2 
y 2 = y" + y'f = a 2 .x i +b 2 .x 2 + a A .x i +b A .x 2 j 
= (t/ 2 + a 4 )Xj + [b 2 +b^)x 2 ' 

Using the matrix notation to represent this same addition 

... ( 2 ) 


(I) 


.»T 

= 

- t . 


'.V’i' 


a x b{ 


x,- 

-h 

'^3 

h a ‘ 

X l 



y" 


y’ii 


a 2 b 2 


x 2 


"4 

/’4 

X 1 


But in the equations (I) above wc have shown that the addition 
of the outputs can be expressed in matrix form as 


(3) 


'JV 


(d, +«.,) (fi, +h 3 ) 1 

X, 

..Vi 


(d 2 +« 4 ) (/>2 + fi 4 ) J 

Xi 


The matrix relations shown at (2) can be expressed 

Y = Y'+Y" 


from the relations (3). 


= A.X + B.X 
= (A + B)X 
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We have demonstrated that the addition of matrices is carried out 
by the addition of corresponding elements of the matrices in the 
correct order and all matrix addition is defined in this way. 

We have also established for matrix algebra that 


A.X + B.X = (A + B)X 


where the curved brackets have the same meaning as in number 
algebra. 

If the circuit M is adjusted so that instead of additional outputs 
being obtained there is a reduction of outputs y ' a and y" f then 


y t - 


y' 


3 a 

yi 




yy 


= AX-BX 


a, b{ 

-xf 


■«3 

bi 

X,‘ 

a 2 b 2 

. X Z. 


a 4 

K 

X 2_ 


(a,-a 3 ) (b,-fc 3 ) 
V(a 2 -a A ) (b 2 -b 4 ) 


= (A - B)X. 


Subtraction of matrices is carried out by subtracting the corres¬ 
ponding elements of the second matrix from the elements of the 
first matrix in the correct order. 

If A — B = 0 we say that the matrices A and B are equal: the 
elements of matrix B must be equal to the corresponding elements 
of matrix A for this to be true. This statement is the definition of 
equal (2 x 2) matrices. 


Examples / 

1 Express the following sets of linear equations in matrix nota¬ 
tions : 

(i) y, = 2x, + 3x 2 (ii) y, = lOx, + 12x 2 

y 2 = 4x, + 2x 2 y 2 = 25x, + 16x 2 


(iii) x' = 3x + 2y 
y' = 2x + 3y 


(iv) x' = 4x-3y 
y' = 2x-y 


2 Rewrite the following matrix equations as linear equations: 


x' 

(i«) 

V 


•3 r 

x' 

y . 




0 2 

.y. 
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II 

u- 

4 6" 

and B = 

2 

3' 


.5 7 


1 

4 


that C = A + B and D = A-B. 


find matrices C and D such 


4 Find with the values of question 3, B + A and compare with 
A + B. 


'2 3' 

, Q = 

1 2' 

, and R = 

1 2 

, find (P + Q), 

1 6 


3 4 


3 1 



(Q + R), (P + Q) + R, P + (Q + R). 

6 With the values of P and Q in question 5, find also Q + P, 
R + Q and compare with P + Q and Q + R. 

7 Using the matrices P, Q and R of question 5, find ( P—Q ), 
(Q-R), (P-Q)-R,-P-(Q-R), and ( P-Q)+R. 


Multiplication of Matrices 


X \ 


y \ 

y\ 



x 2 

M 

-V2 

y 2 

N 

z 2 







Diagram 3 


Now consider an extension to the circuit of diagram 1, so that the 
outputs from network M are fed into N , i.e. become the inputs 
to another electrical circuit N, different from M, from which the 
outputs are z, and z 2 . Then as on page 5 

y, = a 1 .x l +b 1 .x 2 
y 2 = a 2 .x,+fe 2 .x 2 

i.e. Y = A.X with matrix notation. But the outputs from circuit M 
become the inputs for circuit N, so that in a similar manner we 
can state for circuit N 


Z 1 = C,.y,+d x .y 2 
z 2 = c 2 .y l +d 2 .y 2 
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. For circuit N 


Introduction to Matrices 
i.c. Z = B. Y where B represents the matrix 

c 2 (1 2 

are the result of the operation of the coefficients 


< i ‘i i 


the outputs 


<•. </.• 

on the inputs 

wr 

c 2 d 2 


,»2 


operation of the matrix 

a x b x 

on the inputs 

v i 


a 2 b 2 


• Y 2 


To summarize 
and 


Z = BY 
Y = AX 

Z = BY = BAX 


Z l 


<*! 

dv 


~a x 

/>,- 


■*r 

z 2 


C 2 



a 2 

/j 2. 


A‘i 


Note the order of the operations B operates on A . not vice versa. 
By using number algebra we can find by the method of substitution 
the relation between z x% z 2 and x x% x 2 and hence find a meaning 
for matrix multiplication as follows. 

-i = c x .y x -f d\.y 2 and y, = a x .x, + /> t .x 2 

z 2 = c 2 *y\ d 2 • Yz y 2 = a 2 -x x +h 2 .x 2 

By direct substitution for y x and v 2 in the first set of linear 
equations 

f z, = c x (u x .x x +b x .x 2 ) + d x (a 2 .x x +b 2 .x 2 ) 

\z 2 = c 2 {a x ,.v, +h x ,x 2 )+d 2 (a 2 .x x +b 2 .x 2 ) 
f z j = c j. 0 j. x j *F c j. b j. .\*2 4* d j. a 2 * A-d x • b 2 • x 2 
z 2 = c 2 . a |. x | *+* c 2 . b |. ,\*2 4* (1 2 . a 2 . x j 4* d 2 * b 2 • 2 


{ z, = ((', .a, +(/, ,a 2 )x i +(c { ./), + d x . b 2 )x 2 
z 2 = (c 2 .ti x +d 2 .a 2 )Xi +(c 2 ./>i +tl 2 .h 2 )x 2 


(c,£l, +(/,(/ 2 ) ic l b l -f f/,/> 2 )l r.Y, 
(c 2 f/j +d 2 u 2 ) [c 2 b j -Yil 2 b 2 ) x 2 
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i.c. 


but 


U’, .a x +d x .a 2 ) (c*,./?, +il x ,b 2 ) 

Z = 

(c 2 .a x +d 2 .a 2 ) ic 2 .b x +d 2 .b 2 ) 
Z = B.A.X 





(f, +d t .n 2 ) k - i hi +d t .b 2 y 


•* B. A = 





{c 2 .a l +il 2 .a 2 ) ( c 2 .b , +d 2 .b 2 ) 

C, d{ 

a x b x ' 


'(£•, .a, +d, .a 2 ) (c, ./), +d l .h 2 )' 

C 2 cl 2 

KJ 

to 


(c 2 .cii + d 2 .a 2 ) (c 2 .b | +d 2 .b 2 ) 


X 


The elements of the matrices are just numbers and (c x .a x +d x .a 2 ) 
is just a number also. Thus the multiplication of two (2 x 2) matrices 
results in a single (2 x 2) matrix and the form of the process is 
special to matrix algebra. For consistency matrix multiplication of 
(2x2) matrices is defined in such a way that the process can be 
extended to bigger matrices later. Inspection of the products 
above, shows that the elements of a row of the First matrix B are 
multiplied by the corresponding elements of a column of the second 
matrix A, these products being summed for a given row and column 
to form one element of the new matrix. If the first row of B and the 
first column of A are used, the element formed is in the first row 
and first column of B. A : if the first row of B and the second column 
of A are used, the element so formed is in the first row and second 
column of B.A. The row and the column position of each element 
in B.A is thus uniquely determined by the row of B and the 
column of A used to form that element. 

The process of matrix multiplication is easily demonstrated by 
diagrams and quickly learnt by simple numerical exercises. 


Example 


If B = 

1 2' 

and A = 

2 3 


4 3. 


.3 4 


B.A = 

1 2 

rr 


T] 


4 3 

3 


_4_ 


The 1st row' of B is [ 1 2 

The 2nd row' of B is 


4 3 
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1 2 


4 3 


The 1st column of A is 

The 2nd column of A is 


1 

X 

2 


1 

X 

3 

2 

X 

3 


2 

X 

4 

4 

X 

2 


4 

X 

3 

3 

X 

3 


3 

X 

4 


2 


1 2 


3 


3 


4 


2 


4 3 

3 


3 


4 



(1 x 2) + (2 x 3) (1 x 3)+(2 x 4) 
(4 x 2) + (3 x 3) (4 x 3) + (3 x 4). 
'(2+6) (3 + 8) 

.(8 + 9) (12 + 12). 


8 If 
17 24^ 


If we wish to find the product A. B since A is now the first matrix 
we use the rows of + and the columns of B. 


A.B = 


2 3 


3 4 


2 3 


3 4 


T 


2 3 


2 


4 


3 


1 


3 4 


2 


4 


3 
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2 

X 

1 


2 

X 

2 

_3_ 

X 

A_ 


J_ 

X 

_3_ 








~T" 

X 

T1 


T 

X 

T 


X 

4 


4_ 

X 



(2 x 1) + (3 x 4) (2 x 2) + (3 x 3) 

.(3 x 1) + (4 x 4) (3 x 2)+(4 x 3) 

2+12 4 + 9 ‘ 

3 + 16 6+12 

14 13 

119 18. 

At once we notice that this result is different from the previous 
matrix multiplication of B.A 
B.A is not the same as A.B 

A.B * B.A 

The multiplication of matrices is not commutative in general and 
this is at once a very distinguishing feature of matrix algebra. 

The product B.A is called pre-multiplicati on of + by B. 
The product A.B is called post-multiplication of A by B. 
Matrix multiplication can be regarded in another way which has 
an implication for later work with more complicated matrices. 


Again let A = 

’2 3' 

and B = 

1 2' 


.3 4 


4 3 


In matrix A. the elements can be divided into two rows 

Row I R, (R suffix 1 denotes Row 1) 
Row 2 R 2 (R suffix 2 denotes Row 2) 


2 3 

3 4 


In matrix B. the elements can be divided into two columns 
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’rr 


~2~1 

1± 


3 


Column Column 
1 2 

C, c. 


C, denotes column 1 
C 2 denotes column 2 


A.B 



R, c, r,.c 2 
r 2 .c, r 2 .c 2 


In this arrangement the two suffixes on each term denote the row 
and column uniquely. 


J 3j 

l 


2 3 


2 





4 


3 



14 13 

3 4 


i 


3 4 


~2 



sC 

oc 


4 


3 





Examples If 


I If A = 

l 3' 

and B = 

1 2 


6 . 


.3 4 

2 If P = 

1 2' 

and Q = 

- 1 


.3 1 


2 


, find A.B and B.A. 


y 


. find P.Q and Q.P. 


-4 


3 If A = 


2 3' 

. B = 

1 2' 

, P = 

1 

7 

I 6 


3 4 


3 

1 


find A.B then 


multiply this new matrix by P to find (A.B)P: find B.P and then 
find A(B.P). What do you notice about (A.B)P and A(B.P)'? 


4 With the same values for A , B and P as in question 3, find BA. 
(BA)P. AP and B(A.P). What do you notice about [BA)P and 
BiAPYi 
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5 Using the same values for .4. B and P find A + B. P{A + B). 
PA + PB. B + P. A[B + P). AB + AP. What do you notice about 
P[A + B) and PA+PB. about A[B + P) and AB + AP ? 

6 Why does (AB)P of question 3 dilVer from (BA)P of question 4? 

Scalar Multiplication of a Matrix 


k\i 


kx. 



bs 


ky 2 


Diagram 4 


Returning to the original arrangement of diagram 1. it is obvious 
that if both inputs x, and x 2 are increased or magnified by the 
same number k. then both outputs would be altered in the same 
ratio k. i.e. 


kv x = a x ,kx x + h x kx 2 
and 

ky 2 = a 2 .kx x +b 2 kx 2 


k 

wr 


a x .k + b x .k 

•V| 

k 

«, h{ 



.>'2. 


a 2 . k + b 2 . A, 

V 2. 


a 2 b 2 



A ratio number such as k is called a scalar and the above 
shows that multiplying a matrix by a scalar means multiplying 
each element of the matrix by that same scalar. 

With ordinary numbers we have shown earlier that division by 
a number such as 3 is the same process as multiplying by its 
reciprocal j. Division of a matrix by a scalar number such as 4 
is precisely the same as multiplication by the reciprocal of 4, 
i.e. i 

However, the division of one matrix by another is not defined : the 
process is avoided by using the reciprocal matrix in multiplication. 
The reciprocal matrix is now more correctly known as the 
inverse matrix. This process will be developed later. 
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Examples III 


1 If A = 


2 If 8 = 


3 If P = 


2 3 
1 6 

1 2 

3 4 

2 -1 

3 -3 


, find 2/1, 5-4, —3/4. A/1. 

, find 48, kB, -28. 

, find 2P. nP. -3 P. 


If A = 

'4 T 

and 8 = 

'2 3' 


3 2. 


4 5. 


4 Using matrix of question 1. find A + A and compare with 2/4. 

The process of subtraction on page 8 can now be considered 
using the scalar - 1. 

then A—B can be regarded 


-1T2 
4 5 

-2 -3 
-4 -5 


as A + (— 8) or A + (— I )B. 

A-B = 


4 7 

3 2 

4 7 
3 2 


+ 


2 4 
-1 -3 


If 4 = 


In matrix algebra A.A is defined as A 1 . Similarly /4 J = A 2 .A 
= A. A. A. 


'4 7 

then A.A = 

4 7 

4 r 

= 

'37 42 

.3 2. 


3 2 

3 2 


18 25 


The Unit Matrix and Zero Matrix 


The matrix 


I 0 
0 1 


has special properties. 
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2 3‘ 

1 0 


(2 x l) + (3 x 0) 

(2 x0)+(3x iy 

3 4 

0 1 


(3 x l)+(4 x 0) 

(3 x0) + (4x 1) 


2 3 

3 4 


again 


1 O’ 

2 3 


(1 x2)+(0x3) 

(1 x 3) + (0 x 4)' 

0 1 

3 4 


(0 x 2)+(I x 3) 

(0x 3) + (l x 4) 


'2 3‘ 
3 4 


In both cases the multiplication by 
matrix unaltered. This matrix 


1 0 


has left the original 


0 1 

1 01 is called the unit matrix 


0 1 


because it behaves like the number I in ordinary algebra and is 
denoted by /. 


A. I = IA = A always. 

Similarly it can be demonstrated that the zero or NULL 
matrix 0 = TO 0] behaves like the zero of ordinary numbers, i.e. 

0 0 


A C\ _ f\ A 


If we consider any matrix 

post-multiplication by / both 
matrix, thus 


a b 


pre-multiplication by / and 


c d 

yield the same result—the original 


a b m 

1 0- 


a b- 


-1 0 

a b • 

c d 

0 1 


c d 


0 1 

c d 


It is important to note that this unit multiplication holds for 
/ with all other matrices but there are certain matrices where it 
is possible to find another matrix which will behave in a similar 

17 













































































Introduction to Matrices 


manner, thus let A .= 

- 2 0 

and B - 


1 

IJ 

o 



2 01 

1 01 r 


and 


A.B = 


B.A = 


-2 0 


I 2 


_ i 


i 

o 


’ 2 0 


2 Ot 

I 2 


-2 0 


o 

<N 

1 


then 


= A 


= A 


so for these two special matrices 

AB = B.A = A 

and matrix B has behaved in this special case like the unit matrix /, 
but it would not do so with all other matrices as / does. 

Here follows another peculiarity of matrix algebra if 

ri -21 


P = 


PQ = 


3 -6 


and Q = 


r4 2 
2 1 


then 


1 -2‘ 

■4 2‘ 


(4-4) (2—2) 


0 0 

3 -6 

? 1. 


(12—12) (6-6) 


.0 °. 


= 0 


Here we have the product of two non-zero matrices being the zero 
matrix, but 


Q.P = 


Hence in matrix algebra if P.Q = 0 it does not follow that 
either P or Q is zero, or Q.P has zero value. (Here is another 
example of the non-commutative property of matrix multiplica¬ 
tion P.Q # Q.P.) 


-4 2‘ 

1 -2 


(4 + 6) 

(-8-12)- 


10 -20 

2 1 

.3 -6 


.(2 + 3) 

(-4-6) 


5 -10 


Examples IV 


1 If A = 

2 3 

and B = 

1 2 

and / = 

1 0 


.1 6. 


.3 4 


0 1. 


, verify by 


matrix multiplication that Al = lA = A and Bl = IB - B. 


2 If D = 


'3 4" 

, E = 

a b' 

and K = 

‘2 0‘ 

1 6 


c d 


0 2 


, find DK, KD, 


KE, EK. 

18 


Section II 


3 If M = 

2 0 

and B = 

i 

O' 

and / = 

1 O' 


-1 0 


i 

3 


0 1 


MB. BM, IM. and IB. 

4 Using matrices A. B and / of question 1. find 

A-l, A —31. 2B — I, B—kl. 

5 Using the same matrices, find (A - B)(A 4- B). 

Linear Mappings and Transformations 

Transformations using (2x2) Matrices 


3" 

• P'(x',y') 


P(x,y) 


• 

0 

X 

Diagram 5 


If the position of the point P' (co-ordinates x\ y') is related to 
the position of the point P (co-ordinates x, y) by a set of linear 
algebraic equations thus: 


and 


x' = a x .x + b x .y 
/ = a 2 .x + b 2 .y 


then we say that F is a Linear Mapping of the point P. For any 
position of P there is only one possible position of the point F. 
The equations at (1) can be shown in matrix form 


x ,m 


■fi 

b{ 

*X‘ 

y\ 


.°2 

b 2 . 

y m 


and the matrix 


a 


b i 


represents the operation which transforms 


a 2 b 2 

the point P to the position P\ We shall consider some of the 
important transformations and the matrices which cause them. 

P’ is called the image of P. 
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Introduction to Matrices 


Reflection 


Consider the transformation 

x" 


-1 0‘ 

V 


y. 


0 1. 

.y. 



Apply this transformation to two points P(2, 1) and Q( 4, 2) as 
shown, P becomes 



x' 


-1 0 

2 


’—2 

F 

y. 

— 

0 1. 

1 


1 


and Q becomes 


x" 

_ 

-1 0 

■4 


—4' 

y". 


0 1 

2 


2 


It is obvious from the linear relations that the points P and Q 
have been reflected in the y axis to the new positions P' and 

Q". 


The matrix 


-1 

0 


0 

1 


therefore represents a reflection in the 


y axis. 

From similar considerations, the matrix 


1 0 
0 -1 


represents a 


reflection in the x axis. 

In the diagram, F is the image of P and Q" is the image of Q 
and the transformation is similar to the effect of lateral inversion 
produced by a plane mirror because we have used unit values for 
the numbers a, and b 2 . 
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Stretching 



If the points P(2, 1) and Q( 4, 2) are now 


mapping 


3 0 


subjected to the linear 


then P transforms to 


P' 


and Q transforms to 

Q' 


x' 


■3 O’ 

2 


’6’ 

y. 


.0 3 . 

1 


3 

x" 


■3 O' 

'4' 


12 

y. 


0 

2. 


6 


The line bounded by P and Q has been stretched to the line 
P'Q' and the plane has been stretched outwards also. 

If now we combine the two transformations in a single matrix 
— 3 0] this will cause both a reflection in the y axis and a 

0 3j 

stretch at the same time. 


The 4 Shear ’ Transformation 
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Introduction to Matrices 


Consider the four points P( 1,0), 0(4,0), R( 4,2), 5(1,2) and 


apply the transformation 

1 1 

to each 

V 

= 

i i 

X 


0 l . 


y\ 


0 1 

y m 


x' = X + V 

/ = y 


P(l, 0) transforms to P'(l, 0) 
0(4, 0) transforms to 0(4, 0) 
P(4, 2) transforms to P'(6, 2) 
5(1, 2) transforms to 5'(3, 2) 


1 he original points P, 0, P and 5 defined a rectangle of area 6 sq 
units, their images, P', 0', R\ 5', define a parallelogram of the 
same altitude but suffering a shear to the positive direction. The 
area of the new figure P'Q'R'S' also has the area 6 sq units. Hence 
the matrix fl 11 represents a shear of the plane to the right. 


0 1 

The matrix [2 3 
0 2 


would have produced both a shear and 


a stretch of the plane. This is shown in diagram 9. 



Again consider the points P(l, 0), 0(4, 0), R( 4, 2) and 5(1, 2) and 
apply the transformation [2 31 to each. 


P( 1, 0) transforms to P'(2, 0) 
0(4, 0) transforms to 0 (8, 0) 
P(4, 2) transforms to P’(14, 4) 
5(1, 2) transforms to 5 (8, 4) 
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The original points defined a rectangle of area 6 sq units their 
images now define a parallelogram P 0 R S' of area 24 sq units, 
which shows a magnification in both linear directions and a shear 
of the plane also. The ratio of the two areas and the matrix causing 
it will be dealt with in a later section (page 35). 

Rotation of Axes about the Origin 



The point P has co-ordinates a* and y with reference to the axes 
Ox and Ov as shown. PN is the ordinate y and ON the abscissa a*. 
If the axes are rotated through a positive angle 0, to Ox'. Oy 
then PM = y’ in the new axes and OM = x. 

It can be shown by trigonometry that 

x = cos 0.x + sin O.y 

and 

y = - sin 0.x + cos0.y 

i.e. in matrix form 


x' 


cos 0 + sin 0 

V 

y\ 


— sin 0 4- cos 0 

r 


The matrix 


cos 0 sin 0 


represents a rotation of the axes 


— sin 0 cos 0 


through a positive angle 0. This matrix could also represent a 
rotation of OP clockwise about the origin, essentially it is a rotation. 
If the angle of rotation is increased to 90° then sin 0 is equal to 1 
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Introduction to Matrices 


and cos 0 = 0. the 


becomes 


0 

-1 


matrix representing a rotation through 90° 



In diagram 11 the point P(2. 3) is transformed by the matrix 
0 11 to the point P'(3, —2) that is the line OP has been turned 

—! 0 

clockwise to OP’ such that [ POP' =1 Rt.L The same trans¬ 
formation would have been obtained if the axes had been rotated 
anticlockwise through 1 Rt.L Likewise the matrix [0—1 

1 0 


represents a rotation of OP anticlockwise about the origin 0. 

We have now shown that multiplication by a matrix represents 
a transformation so that if this is followed by a second multiplica¬ 
tion by another matrix, another transformation follows the first. 
Now, generally, the order in which transformations are carried 
out alters the final result. We shall show this geometrically. 


Non-commutativc Transformations 



Section II 


If the line PQ is rotated through an angle 0 to position P'Q' and 
then reflected in the y axis, we end with P’Q " being the final 
position. If the reflection is denoted by matrix B and the rotation 
through angle 0 by matrix A % then the operations can be shown 
algebraically. 


and 



for the rotation 


for reflection 


A*"' 


X* 


= B.A 


y" 


y 


A being pre-multiplied by B 


The full transformation of PQ to the final position P 'Q’’ is represen¬ 
ted by the product B.A in that order. Now we shall consider the 
two operations in an alternative order. 



In this case, if the line PO is first reflected in the v axis to position 
P'Q' and then rotated through an angle 0 we end with P 'Q” being 
in a different position from diagram 12. Using matrices A and B 
as before to denote the operations of rotation and reflection in 
the y axis then the operations can again be shown algebraically: 



for the reflection 


R 
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Introduction to Matrices 
and 



for the rotation 


x” 


X 


= A.B 


y" 


V 


But since we have demonstrated that the final result is now different, 
then B.A is not equal to A.B. Once more we have shown that 
matrix multiplication is not commutative. 

The two processes used in diagrams 12 and 13 can be represented 


by the two matrices A = 

cos 0 

-sin 0 

H = 

-1 O' 


sin 0 

cos 0 


0 1 


then in case I 



-1 0 

cos 0 

—sin O' 

B.A = 

0 1 

sin 0 

cos 0 


—cos 0 

sin 0 



sin 0 

cos 0 



cos 0 * 

-sin O' 

-1 0 

A.B = 

sin 0 

cos 0 

0 1 


— cos ft —sin 0 

— sin 0 cos 0 

and because of the sign changes the results arc different, i.c. B.A 
is not equal to A.B (shown B.A ^ A.B). 

Examples V 

1 A linear mapping is defined by x = lx + 3y\ y = 3x + 5y, give 
the matrix of this mapping, express the transformation in a matrix 
equation. 

2 If x ' = 3.v and y' = y, express this transformation in matrix 
form and state what it represents. 
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3 If 

x" 

= A 

X 

and A = 

-1 0 

, how does this transform 


y. 


.y. 


0 1 . 



the points (2, 3) and (3. 5)? 


4 If 

X 

= 

cos 45° 

— sin 45°‘ 

X 

and OP is the line joining 


.y\ 


sin 45° 

cos 45° 

y 



the origin (0, 0) to the point P( 2, 3), show in a diagram the position 
of O P' the image of OP. 


5 Ifrxi- 


Ly 


1 

2 

2 




rx~i. what is the transformation 


produced? 


6 If 

X 

= 

- i i - 

v/2 

X 

, show in a diagram the position of 


- y ' 


1 I 

L-x/2 

-y 



O P' the image of the line OP, joining (0, 0) to P(2, 3). How does 
this transformation differ from that in question 5? 

7 What is the linear mapping produced by the matrix 

8 How does the linear mapping of 


1 1 
0 1 


differ from 


01 ? 
2 
3 

1 


9 Prove by induction that 


cos 0 

- sin 0 

ft 

cos n 0 

— sin n 0 

sin 0 

cos 0 


sin n 0 

cos n 0 
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Introduction to Matrices 


Trigonometrical Relationships 



Let matrix A = 


cos a —sin a 


then the transformation of point 


sin a cos a 

P to P' is represented by A , similarly if B = [cos/? -sin p 

sin ft cos P 

then B represents the transformation of F to P", therefore the 
two rotations which transform OP to 0"P" are given by B.A, 
but the rotation of OP to its final position 0"P" means a rotation 
through an angle (a + /J), i.e., Tcos(a + /?) — sin(a + /?) 

sin(a + P) cos(a + P) 


cos(a + P) - sin(a + P) 
sin(a + /J) cos(a + /?) 


COS P 

— sin p 

cos a —sin a' 

sin p 

cos P 

sin a cos a 


(cos p cos a — sin p sin a) (— cos p sin a — sin P cos a)’ 

(sin P cos a + cos P sin a) (— sin a sin P 4- cos a cos P) 

Equating elements in these two equal matrices we have 

cos(a + P) = cos a cos — sin a sin P 
sin(a + /f) = sin a cos /? + cos a sin p 

These two results are of great importance and can be established 
in other ways. Most students will meet these formulae in a course 
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of trigonometry, but the matrix method is very straightforward 
compared with the other methods commonly used and has the 
advantage over some, that it is general, i.e. applies to angles 
greater than 90° (if the correct signs are used for sines and cosines). 


So far our consideration of matrices has been developed from 
some operation or transformation which has been represented by a 
2x2 matrix, but having given a practical meaning to a 2 x 2 
matrix such matrices can now be considered as an independent 
system without continual reference to any operation which they 
represent. 


The Inverse Matrix 


Let matrix A = 

■5 r 

and matrix B = 

' 3 -7 


2 3 


-2 5 


. Notice that 


matrix B has been chosen to have the same numbers for elements 
but arranged in a different order w'ith some changes of sign 


= I the unit matrix 



5 

T 

’ 3 

-r 


1 

O' 

A.B = 

2 

3 

-2 

5 

— 

0 

1 


again 



3 —7" 

'5 

7 


1 0 

B.A = 

-2 5. 

2 

3. 

— 

0 1 


= / the unit matrix. 


Multiplication between these two matrices is commutative, i.e. 
A.B = B.A = /. 

In number algebra if a.b = 1 then b = ^ = a~ x and wc say 

that b is the reciprocal of u, shown (a l ). 

Similarly in matrix algebra, if A.B = I = B.A, then B is called 
the inverse matrix of A , and is denoted by A~ l . 

A.A^ 1 = A~ l .A = / 


Now lake P = 

'5 3' 

and Q = 

2 -3 


.3 2. 


-3 5 . 
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5 3 

• 2 

: -3 


1 

0 

3 2 

-3 

5 


0 

1 

2 

— 3 

5 3 


1 

0 

-3 

5 

3 2 


0 

1 


P.Q = Q.P = / 

0 is the inverse of P % i.e. Q is P" 


= / the unit matrix 

= / the unit matrix 


In the matrix P = 


[5 31 
\ 

3 N *> 

L * 5 


the number 5 is called the leading 


element and a line drawn from 5 in the top left-hand corner to 
2 in the bottom right-hand corner is called the leading diagonal. 
A line drawn from the bottom left-hand corner to the top right-hand 
corner is called the secondary diagonal. 


Leading elemeni 


It. 


A 

3' An -2 


Sccondar> diagonal 


Leading diagonal 


5 and 2 are the elements of the leading diagonal. 3 and 3 are the 
elements of the secondary diagonal. In each case so far the inverse 
matrix has been formed by transposing the elements of the leading 
diagonal and changing the signs of the elements in the secondary 
diagonal thus 


a b 

INVERSE 

d 

-b 

c d 


—e 

a 


but in each of our examples (ad-be) has been equal to 1. In matrix 


P already used 


V 

3 2 


(5 x 2)-(3 x 3) = 10-9 = 1. 


Now suppose that in the matrix 


b 

d 


, (ad-be) = D. D not 


being 0. and we use a method based on the above for finding the 
Inverse, then providing D is not equal to 0 
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a b ‘ 

• d —b' 


(ad—be) 

( — ab-\-ab)' 

c d 

—c a 


(cd - cd) 

(ad — he) 


D 0 
0 D 


= D 


1 0 

10 1 


= D.l 


Thus, if we now divide each element in the second matrix by 
this value D we have 


a b 

c d 


d__b_ 

D D 

—c a 

D D 


ad —be —ab+ah 


Pre-multiplication gives 
'a b 


t-- 

D D 


— c a 

~D D 


Therefore 


.c d 


d_ -b 
D D 


D 

D 

cd — cd 

1 . 

1 

D 

D 

ad —be 

bd-bd 

D 

D 

— ac+ac 

+ ad—be 


1 0 

0 1 


= / 


1 0 

0 1 


= / 


D D 

will always be the inversi; matrix of 


a b 
c d 


-c a 

~D D 

where D = {ad—be), providing that D # 0. 


It is more convenient to give 


1 


as 


D 


d -b 
D D 
— c a 
~D D 

The quantity (ad—be) is a pure number and is called the deter¬ 
minant value of the matrix 


d -b 
—c a 


The determinant is shown 


a b 
c d 
a b 
c d 


which is distinguished from 
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the matrix by the use of vertical lines on each side of the group of 
numbers. In some books the sign A is used to denote the value 
of the determinant. 

If the determinant value of a matrix is zero, the inverse of the 
matrix does not exist, the matrix is said to be singular. In the 
next section the subject of determinant values and singular matrices 
will be re-examined from a geometrical point of view to give a 
practical meaning to the determinant, but first two examples of 
finding an inverse matrix will be presented. 


Example (i) 


To find the inverse of matrix C = 


2 3 


2 5. 


The determinant D is (2 x 5)—(2 x 3) = 4 
C 


‘=i 

■ 5 -3 

-2 2 

and CC 1 = 

'2 3' 

.2 5. 

1 

4 

5 -3" 

-2 2. 




4 01 

ri oi 


_ i 
- 4 


c~ l .c = i 


= 4 


5 -3 
L-2 2 

4 0 
0 4 

C.C~ l = C“‘.C = / 


0 4 
2 3 
2 5 


0 1 


= / 


= / 


Example (it) 


Let P = 


1 

3 


2 ' 

4 


determinant of P = 4—6 = —2 



P-'.P = 



-2 

1 


2 

4 
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P.P = 

P.P 


= 

1 2 
3 4 

= P l .P = / 


-2 

O ' 


1 

0 

0 

-2 


0 

1 


= / 


(-i) 


4 -21 
-3 1 


= M) 


-2 01 
0 -2 


1 0 
0 1 


= / 


Use of the Inverse Matrix in Coding Methods 

Suppose the letters of the alphabet are assigned arbitrary numbers 

Table 1 


A 

B 

C 

D 

E 

F 

G 

H 

I 

J 

K 

L 

M 

1 

7 

13 

11 

4 

9 

16 

15 

5 

24 

23 

8 

22 

N 

O 

P 

Q 

R 

s 

T 

U 

V 

W 

X 

Y 

Z 

6 

2 

26 

14 

10 

18 

20 

3 

25 

12 

17 

19 

21 


Then the word rent could be coded as 23 4 6 20 and put 
in a 2 x 2 matrix as 


23 4 

6 20 


Let A = 

'23 4' 

and let P = 

'4 5' 


6 20 


.3 4 

If now A 

. is post-multiplied by P the 


be the coding matrix. 


(2 x 2) matrix Q. 
A.P = Q 


The word Kent is now coded in the matrix Q with the aid of the 
coding matrix P, and could be transmitted. At the receiving end, 
if the inverse of P is known, then 


= A(PP~ l ) 

= AI 

= A the original word in number form. 
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Q = 

The inverse of P is 


'23 

4' 

4 

5‘ 


104 

131 

. 6 

20. 

.3 

4 


84 

110 


4 -5 
.“3 4 

P~ l is the decoding matrix. 


because the determinant is 1. 


Now Q.P 1 = 


104 131 T 4 -5 

84 llOj [ — 3 4 

23 4 

6 20 

On reference to the original table the letters represented by the 
numbers give— Kent. 


Examples VI 

1 Find the determinant of matrix 

2 Find the determinant of matrix 


3 If A = 


3 2 


, find 


3 2 


2 3. 

4 5 
4 —6'. 

2 5 . 

What is this class of matrix called? 



. 3 2 


3 

2 



4 If A = 

'3 T 

2 5 

, find A 1 

. Find A A 1 and A~ l .A. 


5 Find the product 

'2 3' 

.4 5. 

r 5 -3' 

.-4 2. 

. Give the inverse of 

'2 3 

4 5 


6 Using the scheme shown in Table 1, put the word deal 
into a (2x2) matrix and encode using the coding matrix 

r 5 3 . Give the decoding matrix P~ 1 and check. 

3 2. 

7 Decode the matrix 


P = 


112 68 

if the coding matrix is 

5 3' 

116 74. 


.3 2. 
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8 

9 

10 


If matrix A = 

If matrix P = 

Matrix Q = 


7 

.2 

P 

t’ 

_ 2 
3 


4 

3 

9 

s 

-5 

8 


. find A *. 

, find P“‘. 

, find Q '. 


11 Find the images P'. Q of the points P(2.1), Q( 4, 3) from the 


transformation due to the matrix A = 


Then find the 


images P", Q" of the points P\ Q' from the transformation due to 
matrix A~ l . What conclusion do you draw? 


12 What transformation is represented by the Inverse of a non¬ 
singular matrix? 


Geometrical Meaning of the Determinant 



If 0. P. Q. R are the vertices of a unit square 0, (0,0) P(1.0), 


P(0, 1), Q( 1. 1) and the>tjansformation due to the matrix 


is carried out then the square OPQR becomes the quadrilatera 
OP'Q'R' where 0(0,0), P'(2. 3), Q'(3. 6), P'(l,3) are new vertices 
as shown. Simple mensuration will show that the area of the 
new quadrilateral is 3 sq. in. and the determinant ot the 
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matrix is also 3. We started with a unit square for convenience, 
but we could have started with any figure such as a triangle and 
the ratio of the new area to the original area would still be 3 for 
this matrix. Always the ratio in which the area is increased is 
equal to the determinant of the matrix which has been used in 
the transformation. 


The inverse of 

2 1 

•s i 

' 3 - 1 

or 

■ i -f 


3 3 


-3 2 


-i §. 


the points 0. P\ Q\ R' undergo the transformation due to the 


inverse matrix 


1 

-1 


-f 

$ 


they go back to the original positions 


defining a unit square. 

The special case of a matrix which is singular follows. 


6 - 

iy 

/*' 

5 - 


/ 

/ 

4 ■ 


fa 

3 * 


✓ 

/ 

2 * 

/p' 

? 

/ 

\ Q 


/ ] 

ip . 

0 

0 1 

2 3 4 5 jc 


Diagram 16 


The matrix 


‘2 1 
4 2 


is a singular matrix (its determinant is 0 and 


it will transform the vertices O, P, Q, R of a unit square to points 
O', P', Q\ R' given by 


x' = 2x+y 


y' = 4x + 2y 

The point O transforms to (0.0) 
The point P(l, 0) transforms to (2. 4)P' 
The point Q(l, 1) transforms to (3, 6)Q' 
The point P(0, 1) transforms to (1, 2)R' 
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The unit square has been transformed into a straight line; the 
quadrilateral OP'Q'R' of figure 15 has now become a straight line 
(in a sense the singular matrix has produced a quadrilateral of 
zero area) but the important result is that the original points 
0. P. Q. R have been transformed to a straight line. This always 
happens with a 2 x 2 matrix whose determinant is zero. The idea 
will arise later with more advanced matrices. 


Now the linear mapping represented by f 2 


4 2 


is not reversible, 


this can be shown by trying to solve the equations of the mapping 

x' = 2x+y 
y' = 4x + 2y 

to get x and y in terms of x' and y we have 

2 x+y = x' 

4x + 2 y = y' 

Attempts to eliminate x or y from the left-hand side fail because 
the determinant 2 1 is zero. 

4 2 

Whenever the linear mapping is represented by a singular 
matrix such a mapping is not reversible. 


In this singular matrix 
elements in the same ratio. 


2 1 
4 2 


the rows and columns have their 


The Singular Matrix 

Now consider the general (2 x 2) matrix 

minant be zero, i.e. ad —be = 0. 

If 

ad —be = 0 
ad = be 
a c 
b d 


a b 
c d 


and let its deter- 
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Introduction to Matrices 


a c 

Suppose the ratio r = k = - then a = kh and c = kd. The 
b d 


singular matrix can now be written 


. For any point P(x % y) 


kh h 
kd d 

x' = kbx + by = b{kx+y) 

y = kdx + dv = d{kx + y) 

For the transform P' of P % 

x b(kx+y) b 

7 ? = ,,, ", , = i (constant) 

y d(kx + y) d 

-♦ x and y' are always in the same ratio and therefore P ' will always 
lie on a straight line. 


We have shown that if matrix 


a b 
c d 


which means that the two rows [a b] and 


a c 

is singular then - = - 
b d 

c d] have their elements 

. , . a c a b 

in the same ratio and if r = - then - = - and the elements of the 
b d c d 

are also in the same ratio. 


two columns 

a' 

and 

h 


c_ 


d 

It follows t 

lat 

if the e 


; ements of the rows, or the elements of 
the columns, of a matrix are in the same ratio, then the determinant 
of that matrix is zero. We can show that the product of a singular 
matrix with any other matrix must always result in a singular 
matrix. 

be singular: we have shown that 


For example, let A = 


this can be written as 


Then if 


a h 
c d 
kb b 
kd d 



v s 


A.B = 

kb b' 

> Q 


(kbp -f bv) (kbq + bs) 


kd d. 

v s. 


.(kdp + dv) (kdq-\- ds ). 
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b(kp + v) b(kq + s ) 
(Hkp + v) d(kq + s) 


Since the elements of the rows or columns are in the same ratio 
therefore C is a singular matrix. Similarly B.A produces another 
singular matrix. 

The inverse of a singular matrix is not defined but we arc 
now in a position to examine another use of non-singular matrices 
in the solution of simultaneous linear equations. 


Examples VII 

I if x' = 4.v + 5v, express this transformation in matrix form. 


>•' = 2.x + 3.v 

Find the images of the points 0(0.0). P(I.O). 0(1. 1) and R(0. 1). 
If O'P'Q'R' is the transformed figure, find its area and compare 
with the unit square. What is the determinant of the matrix of 
the transformation? State the relation between the area of the new 
figure and the determinant. 

2 Sketch the quadrilateral A B CD into which the square ,4( 1.0). 
B(2.0). C(2, 1). 0(1. 1) is transformed by the mapping 


x" 


-3 

T 

x" 

9 

y 


2 

3 

y 


How does the ratio of the areas A B CD' to .4BCD compare with 
the determinant of the matrix f3 2 


3 


Apply the same transformation 


3 2 
2 3 


to the vertices of the 


triangle All. 1). B(4. 1). 0(3.3). Again find the ratio of the area of 
A A B C' with A ABC and compare w ith the determinant of the 
matrix of the transformation. 


4 


Compare the determinant of 


2 3 

3 2 


with the determinant of 


‘3 2" 
2 3 


. Compare the transformations they represent. 
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Introduction to Matrices 


5 In question 1 prove that the ratio of the areas is equal to the 
determinant of the matrix " 


4 5 
2 3 


6 Prove the same result for question 3. 

Solution of a Pair of Simultaneous Equations 


First if we examine the matrix 

minant is (9 — 8) = 1. 

Therefore its inverse matrix is 


3 4 
2 3 


it is seen that its deter- 


3 -4 


and we can apply this 


-2 3. 

to the solution of the pair of simultaneous equations 

3x + 4y = 18 
2x + 3.v = 13 

These can be written 



’3 4 

X 


18- 




2 3 

m y m 


13 



■ 3 -4 

-3 4- 

X 


3 -4 

18 

-2 3 

2 3 



—2 

; 3 

13 


( 54-52) 
(-36 + 39) 


But 


pci rx 

y J [y 

r2 
3 

For these two column matrices to be the same, their elements 
must be equal 

x = 2 
y = 3 
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We would not normally use such a method to deal with a simple 
pair of simultaneous equations but the method can be extended 
to any number of simultaneous equations. This will be dealt with 
later when a general method of inverting any matrix has been 
developed. 

The problem of solving perhaps 100 simultaneous equations in 
100 unknowns could arise in connection with industry, tech¬ 
nological or scientific research. Using classical methods the 
problem could require some months of work by one man. but 
with the aid of a computer and modified method based on matrices, 
such a problem could be solved in seconds. Much technological 
progress has been made possible in recent years by the facilities 
available in computer laboratories. 


Examples VIII 

1 Restate 2x + 3.v = IX in matrix form. 


3x + 5.v = 29 


2 If 

2-3' 

X 

= 

'4' 

, multiply 


.2+3 

X 


6 



the matrices and give the result 


as a pair of simultaneous equations. 


3 if p = 

'2 3' 

, A’ = 

X 

and P - 1 = 

' 5 -3' 


.3 5. 


.y. 


. 

1 


PX, P-'(PX). 


4 Solve 2x + 3y = 18 using the matrix method. 
3x + 5„v = 29 


5 Using the matrix method solve 2x —3y = 4. 

2x + 3y = 6 

6 Solve the equations 3x - 2y = 5 by matrix method. 

8x — 3y = 30 


Z = A.Y and Y=B.X. find 


matrix C where Z = CX. Write down the linear equations which 
Z = A Y. Y = BX and Z = CX represent. 


II 

'2 3' 

, B = 

1 2' 


.3 4. 


.3 -1. 
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Introduction to Matrices 

8 Attempt to solve 3x —2 y = 5 by the matrix method. 

6x-4v = 7 

Explain why the method fails and sketch the graphs of 3x — 2 v = 5 
and 6x — 4y = 7 on the same pair of axes. 

The Transpose Matrix and Orthogonal 
Matrices 

If the rows of a matrix arc interchanged with the columns so that 
the first row becomes the first column and the second row becomes 
the second column then the new matrix so formed is called the 
transpose of the original. 

If A = 


or A' is 


At once it is seen that the leading diagonal is unchanged and 
the determinant of the transpose is also unchanged. 


a h 


c (I 


..then the transpose of A (shown) as A' 


A' = 



11 

1 




b 


d 



Definition 

If the matrix 4 and its transpose 4' have the property that 
A.A' = /. A is said to be orthogonal. 

It can be shown that the only (2 x 2) matrices having this 
property are the matrices of rotation. 

Thus if 



cos 0 

sin 0 


' cos 0 

— sin O' 

A = 

— sin 0 

cos 0 

A' = 

sin 0 

cos 0 


then 

4.4' 


cos 0 

sin O' 

cos 0 

— sin 0 

sin 0 

cos 0 

sin 0 

cos 0 


I 0 
0 1 


/ 


(cos 2 0 + sin 2 (I) 0 
0 (cos 2 0 + sin 2 ()) 


but 


4.4~‘ =/ 
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For an orthogonal matrix 4(2 x 2) the transpose 4' is equal to 
the inverse matrix A -1 . 

If matrix 4' is the transpose of matrix 4, then it follows that 
matrix 4 is the transpose of matrix A'. 

( A'Y = 4 


Definition 


A linear transformation is orthogonal if it preserves the 
length of a vector. 


The rotation represented by the matrix 


cos 0 sin 0 


— sin 0 cost). 

certainly preserves the length of the vector as in diagram 11 and 
is therefore said to be orthogonal. This condition also means that 
its transpose is equal to its inverse. 


Examples IX 


1 If P = 


2 If P = 


2 5 
-2 4 

, find P'. 


. find P' (sometimes written P'). 


3 If Q = 

4 If R = 


2 y 

3 5. 

cos 0 sin O' 

— sin 0 cos 0 
cos x — sin .x 
sin x cos x 
5 To which class of matrices do Q and R belong? 

find 4' and B\ Find (Afl)'. AB 


, find Q 1 and Q'. 
, find R 1 and R'. 


6 4 = 

3 

4 

. B = 

2 3' 


<N 

_i 



r- 


and B'A'. State your conclusion. 

7 In question 6 compare |A| with |A'| and |B| with |£?'|. 

8 Using matrices 4 and B of question 6. show that 


A' + B' = (4 + B)'. 
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Introduction to Matrices 


9 Prove that the only (2 x 2) matrices which are orthogonal are 
the matrices of rotation. 

Complex Numbers 


y 


P'( o,i) 


P" _ 

(-1,0) O 

Diagram 17 


P 

(1*0) * 


It was shown earlier that the matrix 


’ 0 
-1 


1 

0 


represents a clock¬ 


wise rotation through 90°. 

If we consider the line OP when point P is (1, 0) then the opera- 


tion 

0 -1 

T 

gives 

x' 

the new position P' of P, if the 


.1 o. 

0. 


y. 



operation is repeated for this new position then the final position 


of P" will be 


given by 


x"' 


’() -1 

o -r 

T 



1 0 

1 0 

0 


x'" 


-1 O' 

T 


-r 

y" 


0 -1 

0 


0 


The position P"( — 1, 0) shows that the line OP has been rotated 
anticlockwise to OP' since O is transformed into itself. 


But OP is 1 units long and OP' is — 1 units long. 

alters the sign of 


Double transformation 

0 -1 

o -r 


1 0 

1 o 


OP from positive to negative. 
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If we call the matrix 


0 - 1 ' 
1 0 


R, then the operation R 2 is equiva¬ 


lent to the operation - 1. i.e. R is equivalent to V ^-T. But -1 is 
represented in mathematics by i. 


R = 


is the matrix form of i. 


TO -1 
1 0 

This could have been demonstrated differently thus: 


R 2 = R.R = 


O 

1 

1 

O 


- 1 

1 

O 

J o. 

.1 o. 


7 

o 

» 


= (-!)/ 


and again R 2 represents — 1, i.e. R represents N /— 1 or i. 

Leonard Euler first used the sign i (in 1777) to represent the 
square root of minus one in order that he could give the solution 
of the equation x 2 +1 = 0, but the graphical representation of 
y/~-\ as a rotation through 90° was first presented by Caspar 
Wessel in 1797 in a paper on the subject read to the Royal Academy 
of Denmark. Other mathematicians of this period, including Gauss, 

also used the same plan to represent v — 1 graphically and the 
treatment by J. R. Argand in 1806 at Paris was so successful that 
he is remembered most frequently in this work. The argand 
plane to represent complex numbers is a constant reminder of his 
contribution. 

In his work. Gauss realized that the purely imaginary number 
bs/—l and the number which is partly real and partly imaginary, 

a + h v /^T, needed differentiation by name, for the latter kind of 
number he used the word complex which has remained in use 
since. 

Note: Since electrical engineers use the symbol i to represent 

alternating currents they use the symbol j for — 1. The difference 
need not be confusing, in pure mathematics we continue to use / 
but whenever electrical alternating theory is under treatment i is 
replaced by j. 

The imaginary number b s — 1 corresponds to the matrix 


7 

o 


0 

-b 

o 


b 

O 
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Introduction to Matrices 


a 0 
0 a 


, i.e. al , 


The real number a corresponds to the matrix 

so the complex number a+bi corresponds to the matrix given by 
al + b 


0 - 1 


a 0 


o 

1 


= 


+ 


1 0 


0 a 


b 0 


a —b 
b a 


a -b is the matrix corresponding to the complex number 
h a 

a + bi. The leading diagonal gives the real number and the 
secondary diagonal gives the imaginary number. 

The operations carried out with complex numbers correspond 
exactly to the matrix operations w'ith the matrix form. 

Consider addition of a + bi and c + di. Adding the real and 
imaginary parts separately gives 


(a + hi) + (c + di) = (a + c) + (h + d)i. 

The corresponding matrix operation produces a similar result 


a —b' 

+ 

c -d 


'(fl + C)- 

-( b+dy 

b a 


d c 


(b+d) 

(«+(•) 


Again multiplication of the same tw'o numbers corresponds 
exactly: 


a —ft’ 

c -d- 


(ac—bd)—(ad+bcY 

b a 

d c 


(be + ad) (ac — bd) 


Real part ac—hd , imaginary part (ad + bd)i and 
(a + bi)(c + di) = ac + i 2 bd+(bc + ad)i 


= (ac — bd) + (bc + ad)i since i 2 = — 1. 


So the group of 2 x 2 matrices of the form 


a 


-b 


gives exactly 


b a 

corresponding forms for addition and multiplication as the group 
of complex numbers of the form a + bi. This exact correspondence 
between tw'o groups is called isomorphism and the group of 
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(2x2) matrices which correspond to the complex numbers is said 
to be isomorphic to the group of complex numbers. Many 
problems in complex numbers can be solved by the use of the 
corresponding isomorphic forms of matrices. Thus 


Example (t) 


(3+20(2 + 3/) becomes 

'3 —2' 

’2 —3' 



2 3 

3 2 



= 13 

(3 + 20(2+3/) = 13/ 


0 -13 
13 0 

0-1 
1 0 


Example (//) 

To find 3 + 2/ 

2-3/* 

Since division in matrix algebra is not defined w^e must multiply 
by the inverse matrix. 




The matrices \a — /?] where a and b take all real values can be 


considered as the mathematical model of complex numbers. 


‘3 

-1 

13 

3 1 
“13 


0 

131 

“13 


0 

-r 

2 

3 

.13 

2 

T3J 


13 

Ll3 

0 


i 

0 


Examples X 

1 Add 2 + 3/ and 3 + 2/ using matrix form. 

2 Subtract 2+3/ from 6 + 4/ using matrix form. 

3 Find (2 — 3/)(4 — 5/) using matrix form. 

4 Find (2 + 3/) 2 using matrix form. 

= (2 + 3/K3 — 2/)” *, hence, using the matrix forms, find 

2 + 3/ 

3 — 2/* f 

6 Put r cos 0 + r sin Oi in matrix form, r 
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Introduction to Matrices 

7 If Zj = r,(cos 0! + / sin 0 t ) and Z 2 = r 2 (cos 0 2 + / sin fl 2 ), ex¬ 
press each of these in matrix form and then find Z { .Z 2 and simplify. 


Vectors 

On page 5 single column matrices were called column vectors 
or column matrices. We are now' at a stage when this idea can be 
expanded further. 

The engineer and scientist have long recognized quantities which 
have both magnitude and direction and other quantities which have 
only magnitude or size. Quantities having magnitude and direction 
have been given the name vector. The name was adopted by Sir 
W. R. Hamilton about 1853 when he first published a full treat¬ 
ment on Quaternions. Quantities w'hich have magnitude only can 
be represented by a single number and are called scalars. 

The engineer has a very clear idea of w'hat he means to be 
vector quantities, like force, velocity, momentum, 
acceleration. The mathematician finds that he can represent 
vectors in a plane by a pair of ordered numbers, vectors in three 
dimensional space by a triple of ordered numbers and mathe¬ 
matical quantities which require four or more ordered numbers 
to specify them are for him an extension of the system of vectors. 
A mathematical quantity which requires n ordered numbers to 
specify it completely is called a vector of n dimensions. 



If OP represents a vector in a plane, the distance OP represents 
the magnitude and the direction is given by the angle 0 and the 
sense by the arrowhead. The complete direction is given by the 
two ordered, real numbers (x, y) since y/x = tan 0 and the sense 
O to P is determined by the algebraic sign of x, and y. If P' is 
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the point (-x, y) then the new' sense is shown, the magnitude 
remaining constant. By Pythagoras the length OP = y /x 2 +y 2 
shown mod. OP or \OP\. 


Our earlier use of 


was in fact then a vector representation 


in a space of two dimensions, hence the use of the term column 


vector. Now the transpose of 


is defined as [x y] called a 


ROW VECTOR. 

On page 19 w'here a linear mapping is explained the trans¬ 


formation which operates on the column vector 


y 


produces 


another column vector 
matrix form 


so that in this case the equation in 


.. ( 1 ) 


V 


«i 

V 


X 

y. 


Sh 

K 


y. 


represents a vector relation. 


The matrix 


a, 

a 2 b 2 


can be regarded as the composition of 


two column vectors or of two row vectors. 

In the relation (1) above, the row vector or row matrix [a, 6,] 


operating on the column matrix 


gives the single equation 



II 

since 

x' = a,x + b,y 

and 

y'= [a 2 M' 

i.e. 

y = a 2 x + b 2 y 
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Introduction to Matrices 



Let the vector OP be specified by the ordered pair of numbers 

(x, y). Then the vector (Ax. Ay) will have the direction of OP. but 
depending on A, OP' will fall along OP and if A is greater than I 
there will be magnification. If A is less than 1, there will be com¬ 
pression but mathematically this is still termed magnification. 

If P' is the point Ax, Ay then vectors OP and OP' have the same 
direction (or are parallel). 

If A is a negative number the vector OP' will have the same 
direction but opposite sense to OP. 

The engineer working with vector quantities regards the vector 

OP as the resultant of two displacements ON to the right and NP 
upwards as shown. All the properties of vectors which the engineer 
and scientist use are included in the mathematical model which 
the mathematician has built with the system of ordered pairs of 
numbers, or ordered triples of numbers. 


The Addition of two Vectors 



We use the column matrix 


to represent a vector. In the 
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diagram OP' and OP” are two vectors OP’ = 


x" 

and OP" = 

x" “ 

/. 


y" 


OP' + OP" = 
becomes 


by the method of matrix addition this 


x' + x" 

y'+y" 


In the figure, OL = ON + NL = ON + P'K = x' + x" 
AL = AK + KL = y"+y’ = y'+y" 


The vector OA = 


x'+x" 

y'+y" 


and represents the addition of OP' 


and OP". 

This form of addition of two vectors agrees exactly with the 
form adopted by the engineer and scientist when he is adding 
vector quantities. 

Examples XI 

—> —► 

1 (i) Show diagrammatically the vectors OP and OP' where 

OP = 


(ii) OP + OP = 
matrices = 


3" 

and OP' = 

r 

2 


4 


= 

T 

+ 

r 


2. 


4 


= (3+1). By addition of column 
(2 + 4) 


= OR. 


On your diagram show OR the sum of the two vectors 
OP and OF. 

2 Using the same vectors OP = f3l OP' = \ 1 show 


3 

OP' = 

1 

.2. 


.4. 


OP-OP' = OS, 

where OS is obtained first by subtraction of column matrices. 
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Introduction to Matrices 
3 Show the vectors OA = 


4 

OB = 

0 

0 


3 


and OC = OA + OB s 


give vector OC as a column matrix. 


4 


Show the vector OP = 

'4 

OQ = 

2' 


2 


5 


and OR = OP + OQ. 


Join the points PR and QR. Which kind of quadrilateral is formed 
by OPRQ1 


5 Draw a pair of straight lines OA and OB as axes at any acute 
angle AOB (say 70°). Now show the vectors of question 4. What 
is the figure OPRQ in this case with oblique axes? 


Matrix as a Store of Information 

The section on coding methods on page 33 shows that a matrix 
can also serve as a store of information classified in numerical 
form. This illustrates one of the most important uses of matrix 
methods in linear Programming, Computer Control and in the 
analysis of many processes such as Progranvned Instruction Studies 
and Psychological Testing . Four letter words can be stored in a 
2x2 matrix by the method shown: processes could be stored in 
matrix form by a similar coding. Other matrices of a special form 
can then operate on the stored information. Some of these special 
matrices are given. 


Matrices of special operation 


The matrix 


0 1 
1 0 


causes the operation of reversal thus 


(i) Post-multiplication causes interchange of columns: 


a b' 

0 

1- 


'/> a' 

c d 

1 

0 


d c 


(ii) Pre-multiplication gives interchange of row: 


0 

1- 

a 

b- 


c 

d- 

1 

0 

c 

d 


a 

b 
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causes the operation of addition thus 


(i) Post-multiplication 
(addition of columns) 

(ii) Pre-multiplication 
(addition of rows) 


a 

b' 

l 

r 


a a + b' 


c 

d 

0 

l 


c c + d 


1 

r 

a 

b‘ 


a + c h+d 

0 

i 

c 

d 


c i 

d 


The matrix A' = 


1 01 
1 1 


is the transpose of the matrix 


A = 


1 1 
0 1 


i.e. rows and columns have been interchanged. 


The transpose 


1 0 
1 1 


also causes addition thus 


(i) Post-multiplication 

a b~ 

-i 0 


a + b b ' 


c d 

1 1 


c+d d 


But the addition is now stored in the first column. 


(ii) Pre-multiplication 

1 O' 

a b 


a h 

shows that 


1 1 

c d 


a+c h+d 



the addition is now stored in the second row. 


The two matrices 


i 

1 

and 

i 

0 

0 

1 


i 

1 


are both examples of 


triangular matrices having the properties illustrated above. 
They are called triangular because the elements lie to one side 
only of one of the diagonals. 


ri ii 

and 

1 01 

\ 


\ 

o \ 

f 

.1 \ 




These forms will be met again in more complicated forms of 
matrices. 
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Introduction to Matrices 


Elementary Matrices 

The types of operations on rows (or columns) such as the addition 
of rows or reversal of columns or interchange of rows are pro¬ 
duced by elementary matrices. 

The five elementary matrices and their geometrical equivalent 
interpretation: 


I [0 
.1 0. 


A reflection of the plane in the line through 
the origin at 45° to the x axis. 

ii r« o' 

.0 1 


An elongation or compression parallel to the 
.v axis. 

m ri o' 

.0 a. 


An elongation or compression parallel to the 
y axis. 

iv ri oi 

b 1_ 


A shear parallel to the y axis. 

v ri b' 

o 1. 


A shear parallel to the x axis. 


Any linear transformation of the plane may be represented as 
a product of reflections, shears and elongations in one dimension. 


Nilpotent Matrices 


A matrix of the form 


0 2 
0 0 


has special properties also. 


If 7 = 

0 2 

then 7.7 = 

0 2' 

0 2' 

= 

0 0 

• 

0 0 


.0 0. 

.0 0. 


.0 0. 


. the null 


matrix 


Le. 7 2 = 0 

=> T 3 = T 2 .T = O.T = 0 
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Thus T is a matrix which has the special property namely that 
powers of itself are zero. 

In number algebra the only number whose square is zero is 
zero itself. 

Here then is another difTerencc between matrix algebra and 
number algebra. The property of matrices like T will be met again. 

Such matrices as 7 are said to be nilpotent. 


Examples XII 


I A = 


2 B = 


3 C = 


4 D = 


5 7 = 


6 7 = 


3 4 
-2 I 

4 3 
2 - 2 . 

5 -2 
1 3 
10 2 

3 -6 


. H = [0 I 
.1 0 
and K = 


and L = 


. find AH and HA. 

. find KB and BK. 

, find LC and CL. 


1 3 
0 1 
1 0 
-2 1 


and N = 

' 1 0 


.0 3. 


find DN and ND. 


0 3 
L0 0J 
0 2 
0 2 


. A = 3 4 

.-2 I 
and B = 


4 3 

•) _•) 


II 

a h 

. A = 

-2 O' 

. B = 

«*■ i 

O 


.c cl 


0 -2 


• 

o 


find T 2 and TA and AT. 

. find F 2 . FB and BF. 

. find AP. PA. AB. 


BA. BP. PB. 

8 What geometrical transformations are represented by the 
matrices H. K. L and N above? 


Distributive Laws for Scalar Multiplication 

Scalar multiplication of matrices will now be extended. We have 
demonstrated earlier that the multiplication of a matrix by a 
scalar number k is equivalent to multiplication of each element of 
the matrix by k. 
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Introduction to Matrices 


Let A = 




and 

B = 

b x 

V 



.«3 

« 4 . 





h 4 . 


kA = 

ka 1 

ka 2 

and 

kB 

= 

A 7 >, 

kb 2 


ka 3 

ka 4. 




,kh 3 

kb A . 


kA + kB — 

ka x ka 2 

+ 

kb t kb 2 " 


. ka 3 ka 4 . 


kb 3 kb A 


ka x +kb x ka 2 + kh 2 
. ka 3 + kh 3 ka A + kb 4 _ 
lda x +b x ) k(a 2 + b 2 )' 
„k(a 3 + b 3 ) kUi^ + hi). 

\(a x +b x ) (a 2 + b 2 ) 

= k 

,(a 3 + b 3 ) (04 + ^)- 
= kiA + B) 


Similarly it can be shown that k(A + B) = kA+kB. These two 
relations are the Distributive Laws for Scalar Multiplication. 


Laws for the Addition and Multiplication of 
(2 x 2) Matrices 

At this stage it is useful to summarize some of the work done so 
far. At the beginning we showed by practical situations involving 
operations or transformations that (2x2) matrices could represent 
such processes. It later became obvious that (2 x 2) matrices 
constitute an independent system with an algebra of its own. 
similar in many ways to ordinary number algebra, but different 
in one very important feature that multiplication is not com¬ 
mutative. Such a system of (2x2) matrices constitutes a ring. 
From the examples given the algebra of (2 x 2) matrices obeys the 
following laws for addition: 

(i) A + B = B + A (Commutative Law for addition) 

(ii) (/1 + B) + C = X-HB + C) (Associative Law for addition) 

(iii) A + O = A (Null or Zero matrix is the Identity element for 
addition) 
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where the + sign has the special meaning in matrix algebra of 
adding corresponding elements. 

The multiplication of (2 x 2) matrices follows these laws: 

(i) ( AB)C = A{BC) (Associative Law for multiplication) 

(ii) M + B)C = AC + BC (Distributive Law for multiplication) 
C(A + B) = CA + CB Note the order of multiplication since 

multiplication is not commutative 

(iii) AI = IA = A Unit matrix is the Identity element for 

multiplication. 

Where multiplication, shown A.B or simply AB , has the special 
meaning of matrix multiplication, it must always be emphasized 
that although we use the familiar symbols of number algebra, 
nevertheless matrix algebra is a new algebra with new meanings 
given to the old symbols. In matrix multiplication the position of an 
element in row and column is of the utmost importance and for 
many purposes it is desirable to use a system to identify the row 
and column of an element. This can be done by adopting the con¬ 
vention that two suffix numbers are used so that the first suffix 
number identifies the row and the second suffix number identifies 
the column. 

w, 2 "1 un means the element in row 1 and column I 
a 21 a 12 1 a i\ means the element in row 2 and column 1 

The suffix number system has great advantages particularly 
in later work, but for a beginner the method used so far has other 
advantages. 

When we come to regard matrices as a system of independent 
operators which can be regarded as having their own algebra 
without reference to the operations they represent, then the general 
definition of a matrix as an array of numbers is obviously accep¬ 
table. This definition was due to syi.vf.sthr a contemporary of 
cayley and these two men arc always remembered as the founders 
of modern matrix algebra and matrix notation. 

In the next section we shall extend the ideas already developed 
to larger matrices that deal with mathematical quantities which 
require three ordered numbers to specify them, i.e. vectors of three 
dimensions. 


c 
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Section III 

Meaning of a (3 x 3) Matrix 


X| 

H 

y 1 

X 2 

yi 

*3 

y 3 





Diagram 21 


Consider the network represented by H with three electrical inputs 
Xjl, x 2 , x 3 (these inputs could be e.m.f.s. or currents or a combination 
of both). The total input to H is a mathematical quantity specified 
by the three numbers x u x 2 , x 3 * i- e - the input is a vector of three 


dimensions shown 


*2 

L* 3 J 


as a column matrix. If the outputs 


y t , are eac h a function of the input 


relation is linear then 


*2 

lX 3 j 


such that the 


y i = a l x l + b l x 2 + c l x 2 


y 2 = a 2 x x +h 2 x 2 + c 2 x 2 


y 3 = tf 3 X l+M 2 + ^ 3*3 


i.e. 


> r 


a, b t c,' 


*r 

y 3 

= 

a 2 h 2 c 2 


X 2 

.J 3 - 


0 3 h 3 c 3 . 


_ X 3_ 
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If Y is the column matrix 

~y r 

X is the column matrix 

•V 


i'2 


X 2 


.yi. 


-* 3 _ 


A now represents the matrix above, which has three rows and three 
columns, then 

Y = A.X 

A is called a (3 x 3) matrix and the column matrix X multiplied 
in matrix manner by the (3 x 3) matrix A gives the column matrix Y. 
Again the matrix represents the operation of the electrical net¬ 
work H on the inputs x 1# x 2 , x 3 converting them into the outputs 

yi*y*y* 

Addition of (3 x 3) Matrices 


X \ 

H 

y 

*2 

>•" 

*3 

y'" 




Diagram 22 


Now suppose that the network H is adjusted to give a small output 
/, y"\ again depending on the inputs according to the linear 

relations 

y' = a,x 1 +h t x 2 + c 1 x 3 
y" = a 2 x l +b 2 x 1 + c I x 3 
y‘" = a 3 x t +b 3 x 2 + c 3 x 3 
In matrix form this is 


-y" 


'fli 

K 



y" 

= 


b 2 

<2 

*2 

_y’"_ 


«3 

b 3 

‘• 3 - 



or T = AX. 


59 















































Introduction to Matrices 

Next suppose that an additional output is obtained calling these 
additional outputs y' v yj, y"\ we then have 

= d x x t +c,x 2 +/,x 3 

y'fi = d 2 x x +e 2 x 2 +f 2 x 3 

y"! = d 3 x t +e 3 x 2 +f 3 x 3 


or 


Y" = 


d x e x J\ 
d 2 h 

di f 3 J 


= B.X 


If the two outputs 

’ y' ~ 

and 

■ 


y" 


y'h 

or total output f y x 

_y'"_ 

y. 



X\ 

• Y 2 

L*3J 

are now added to give the final 


y2 

U’3 

y , = y'+y = (a t x i + b i x 2 + c l x 3 )+(d l x { +e l x 2 +f l x 3 ) 
>’2 = y"+y't = ( fl 2 x i + b 2 x 2 +c 2 x 3 )+(d 2 x t +e 2 x 2 +f 2 x 3 ) 
y 3 = = (a 3 x t + b 3 x 2 +c 3 x 3 )+(d 3 xi+e 3 x 2 +f 3 x 3 ) 


y i 
>’2 
y 3 . 


= (a, +d,)x, +( 6 , +e,)x 2 +(c, +/,)x 3 
(a 2 + c/ 2 )x, +( 6 2 + e 2 )x 2 +(r 2 +/ 2 )x 3 
(n 3 + d 3 )x t +(b 3 + e 3 )x 2 +(c 3 +/ 3 )x 3 


y = 


(a,+</,) (6,+e,) (e,+/,)' 
(n 2 + </ 2 ) ( 6 2 + c 2 ) (c 2 +/ 2 ) 
(a 3 + </ 3 ) ( 6 3 + e 3 ) (c 3 +/ 3 )_ 


(I) 


But y = r+ y" = AX + BX = (/1 + B).Y from relation (I). There¬ 
fore as with ( 2 x 2 ) matrices so with (3x3) matrices, addition of 
two (3x3) matrices is carried out by the addition of corresponding 
elements and further the process shows again that 

AX + B.X = (A + B) X 
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Scalar Multiplication 

If B = A. i.e. the elements of B are equal to the corresponding 
elements of A then it follows that 


A + A = 

'2 «, 

26 , 2 c ,1 



2a 2 

: 26, 2 c 2 



2a 2 

, 2b 3 2c 3 . 


=> 2A = 

* 2a , 

26, 2 c,* 



2 a 2 2b 2 2c 2 



.2a 3 2b 3 2 c 3 _ 


Similarly 3A = 2A + A 

= 

‘3^! 3 b x 3 c |" 



3a 2 3^2 3c 2 



_ 3fl 3 3/? 3 3c 3_ 

The value of k.A = 

’ ka 

, kh { kc x 

• 


ka 

2 kb 2 kc 2 



_ka 

3 kh 3 kc i. 



k is a pure number in the system and is called a scalar; kA is 
called scalar multiplication. 

Subtraction of (3 x 3) Matrices 

If in the network H the first output Y' had been followed by 
a reduction of output Y" then the final output would have been y. 

y = Y'-Y" 


= AX — BX 
= (A-B)X 


and the process of subtraction would have given also 
'(a,-r/,) ( 6 ,-c,) (o,-/,)* 
y = (a 2 —d 2 ) ( h 2 -c 2 ) (c 2 -f 2 ) X 
L (^3 — d 3 ) ( b 3 — c 3 ) (c 3 fi ,) 
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Hence subtraction of two matrices is carried out by the subtrac¬ 
tion of elements of the second matrix from the corresponding 
elements of the first matrix. 

=» A — B = A+( — B) 

= A +(- I B) where — 1 is a scalar. 

Subtraction has been shown as a process of addition of negative 
elements. 


Examples XIII 

1 Express in matrix form 

= 2x l + 3x 2 + 4x 3 
y 2 = x,-2x 2 + 3x 3 
)'s = 3x,+2x 2 —2x 3 


2 Express as linear equations the matrix equation Y = AX, 


where A = 

' 3 

2 

r 

. Y = 

"y i" 

. x = 

**i" 


-2 

3 

6 


>’2 


x 2 


I 

-2 

4. 


J’3- 


-■*3- 


3 If P = 


Q-P. 

4 If B = 


5 If P = 


1 2 1 

3 2 1 

2 3 2. 

1 2 r 

3 2 3 

2 I 2 


and Q = 


3 2 

4 2 


1 3 1 


, find P+Q, Q + P, 


find 3 B and kB. 


'i 2 r 

II 

Oi 

i i i" 

. R = 

I 

o 

1 0 1 


2 1 2 


2 1 3 

.3 1 2. 


.3 2 3. 


.1 2 3. 


find 


P + Q, (P + Q)+R, ( Q+R), P+(Q + R). What can you say about 
(P + Q) + R and P+(Q + R )? 
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Multiplication of (3x3) Matrices 





y i 


fl 

*2 

Af 

^2 

3^2 

N 



3^3 

>3 

f3 







Diagram 23 


Consider two separate circuits M and N so that the output of M 
becomes the input to N. If the operation of circuit M transforms 
the inputs x„ x 2 , x 3 into the outputs y„ y 2 , y 3 which in turn are 
transformed by circuit M into final outputs Zj, r 2 , z 3 , then using 


the notation of diagram 21 

• 

y, = a 1 x,+b 1 x 2 + c,x 3 

z, = d x y x + e x y 2 +f x y 3 

y 2 = a 2 x t + b 2 x 2 + c 2 x 3 and 

z 2 = d 2 y\ A~ e 2 y 2 +f 2 y 2 

y 3 = a 3 x,+f> 3 x 2 + c 3 x 3 

23 = ^3^1+^3> , 2+/3>’3 

In matrix notation 


Y = A.X and 

Z = B.Y 

- Z = B(AX) = 

B.A.X 


but in algebraic notation 

Z = z, = ^(fltX^fcjXj+CiX^+etfajX^ftjXj+CjXj) 
+fi(a 3 x l + b 3 x 2 +c 3 x 3 ) 


z 2 = d 2 (a l x l +b l x 2 +c l x 3 )+e 2 (a 2 x l +b 2 x 2 +c 2 x 3 ) 
J 2 [(t 3 x | + b 3 x 2 + C 3 x 2 ) 


z 3 = d 3 (a i x l +b l x 2 + c i x 3 ) + e 3 (a 2 x l +h 2 x 1 + c 2 x 3 ) 
+f 3 (a 3 x 2 +b 3 x 2 + c 3 x 3 ) 
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-2 


z = T',1 = ld x a l +e l a 2 +f l a 3 )x i +(d l b l +e x b 2 +f l b 3 )x 2 

+(</,£•,+?, C 2 +f l C 3 )X 3 

= (d 2 a x + e 2 a 2 +f 2 a 3 )x x +(d 2 b , + e 2 b 2 +f 2 b 3 )x 2 
+(</ 2 c, +e 2 c 2 +f 2 c 3 )x 3 

= (d 3 u x +e 3 a 2 +f 3 a 3 )x x + (d 3 b x +e 3 b 2 +f 3 b 3 )x 2 
+(d 3 c ] +e 3 c 2 +f 3 c 3 )x 3 

adopting the scheme shown on page 14 this gives 


dt e x f x 


d i e x /, 




b 3 


d i e x /, 


di e 2 ii 


d 3 e 3 


M 


a 3 


d 3 e 2 J2 


d2 e 2 fz 


d 3 e 3 f 3 


d 3 e 3 f 3 


Xz 


L X 3 J 


(Row, of BXCol, of A) (R, 

ofBXC 2 of/1) (R, 

ofPXC 3 of/4)" 

"*l‘ 

R 2 C, 

R 2 C 2 

R 2 c 3 

-\'2 

R 3 C 1 

r 3 c 2 

R 3 c 3 

-^ 3 . 


Hence multiplication of (3x3) matrices follows exactly the 
same form as matrix multiplication of ( 2 x 2 ) matrices, this being 
illustrated by a numerical example: — 


r l 

2 

1 ‘ 

Q = 

"3 

2 

f 

3 

2 

1 


1 

1 

2 

.2 

3 

2 . 


.1 

2 

1 . 
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If we multiply in the order P.O then we consider rows of P and 
columns of Q as with (2 x 2) matrices: 


P = 




R,, R 2 . R 3 rows of P C,. C 2 . C 3 columns of Q 


then from the relation just established. 


P.Q = 

R,C, 

r,c 2 

R ,c 3 


R 2 C, 

R 2 C 2 

R 2 c 3 


.RjC, 

R 3 c 2 

R 3 C 3 . 








— 

1 2 1 


3 


1 2 1 


2 


1 2 1 


1 



1 


1 


2 



1 


2 


1 


3 2 1 


3 


3 2 1 


2 


3 2 1 


1 



1 


1 


2 



1 


2 


1 


2 3 2 


3 


2 3 2 


2 


2 3 2 


1 



1 


1 


2 



1 


2 


1 
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Introduction to Matrices 


1x3 

2 x 1 
I x 1 
3x3 
2x1 

1 x 1 

2x3 

3 x I 

2 x 1 


'3+2+1 

9+2+1 


L 6 + 3 + 2 

' 6 

6 

12 

10 

.11 

11 


1x2 
2x1 
1x2 
3x2 
2x1 
1x2 
2x2 
3 x 1 
2x2 


2 + 2 + 2 
6 + 2 + 2 
4 + 3+4 

6' 

8 

10 . 


1 x 1 
2x2 
I x 1 
3x1 
2x2 
1 x 1 
2x1 
3x2 
2x 1 


i+4+r 
3 + 4+1 
2 + 6+2. 


Notice that the multiplication of a (3x3) matrix by another 
(3 x 3) matrix results in a (3 x 3) matrix also. If we wish to repeat 
this process of matrix multiplication but post-multiplying Q by P. 
i.e. to find Q. P, then the elements of Q must now be arranged in 
rows and the elements of P will in this case be arranged in columns 
as shown 
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R,C, 

r,c 2 

R.c 

R 2 C, 

r 2 c 2 

R,C 

.RaC, 

RjC 2 

R 3 C 


We can now proceed to the step showing the multiplication of 
elements. 

3x2 
2x2 
1x3 
I x 2 
1x2 
2x3 
1x2 
2x2 
1 x 3 


3x1 

2x3 

1x2 

1 x 1 
1 x 3 
2x2 
I x 1 
2x3 
1x2 


3x1 

2x1 

1x2 



l x 1 

1 x 1 

2x2 



1 x 1 

2 x 1 

1x2 



3+6+2 6+4+3 
1+3+4 2+2+6 
1+6+2 2+4+3 


3 + 2 + 2' 
1 + 1+4 
1 + 2 + 2 . 


11 13 7 

8 10 6 

9 9 5. 


It is obvious that the resulting (3x3) matrix from Q.P is not the 
same as P.Q. 

Once more the non-com mutative property of matrix 
algebra is illustrated for multiplication, but addition does follow 
the commutative law and we can show this with the two matrices 
P and Q. 

*1+3 2 + 2 1 + 1' 

3+1 2+1 1+2 
. 2+1 3 + 2 2 + 1 . 


P + Q = 


67 



















































































Introduction to Matrices 


'4 4 2" 

4 3 3 
-3 5 3. 

But since the elements are numbers the order of addition is not 
important, hence 


Q + P 


and therefore P + Q = Q + P. 


'4 4 2 
4 3 3 
.3 5 3. 


The Null or Zero Matrix 

As with (2 x 2) matrices the null or zero matrix for (3 x 3) 


matrices is 


0 0 0 
0 0 0 


= 0. 


.0 0 0 . 

Addition of elements gives P + 0 = P 

and 0 + P = P 

Multiplication of P by 0 must result in zero elements and hence 
the product is a zero or null matrix 

P.0 = 0.P = 0 


The Unit (3 x 3) Matrix 

With (2x2) matrices, the array 


1 0 
0 I 


behaves as the unit 


matrix, it is reasonable to suppose that the matrix 


behaves in the same way with (3 x 3) matrices. 
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1 0 0 
0 I 0 
0 0 1 


Section III 


This can be tested with a numerical matrix. 


Let A = 


1 

2 

I 


2 3 
I 1 

3 2 


12 3 
2 1 1 
1 3 2J 

Post-multiplication by 


1 0 0 
0 1 0 
0 0 1 

1 0 0 
0 1 0 
0 0 1 


1 2 3* 

2 11 = A 
1 3 2. 

has left A unchanged. Again, 


pre-multiplication gives 


The matrix 


'1 

0 

o' 

‘1 

2 

3' 


"i 

2 

3' 

0 

1 

0 

2 

1 

1 

= 

2 

1 

1 

.0 

0 

1 . 

.1 

3 

2 


.1 

3 

2. 

ix 

1 

0 

O' 

has 

behaved 

like 

the 


0 

1 

0 








.0 

0 

1 . 








= A 


number algebra, the matrix is called the unit matrix for (3x3) 
matrices and again is denoted by I. 

A.I = 7/4 = 4 


Taking the general matrix M = 


a, b, c | 
a 2 b 2 c 2 


«3 h s 


■3J 


post- or pre¬ 


multiplication by / always yields an unchanged matrix. 


"“I 

b , 

c, 

~1 

0 

O' 


’«1 

b\ c,' 

a 2 

h 2 

C 2 

0 

i 

0 

= 


b 2 c 2 

-“3 

l> 3 

Cl. 

0 

0 

1. 


.«3 

b 2 c 3 _ 


= M 
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"1 

0 

O' 


b x 

<T 


a \ 

b x 

C t ~ 

0 

1 

0 


b 2 

<■2 

= 

“2 

b 2 

<*2 

.0 

0 

1 . 

-«3 

b 3 




b 3 



M . / = / . M = M for all the different Ms. 


The Scalar Matrix 


The matrix 


k 0 0 
0 k 0 
0 0 kj 


is called the scalar matrix and is 


obviously obtained by the scalar multiplication of matrix l by k 


kl 


~k 0 0" 
0 k 0 
0 0 k 


The Diagonal Matrix 

The scalar matrix is a special case of a diagonal matrix. 
In a diagonal matrix all the elements except those on the leading 
diagonal are zeros. 


2 0 0 
0 3 0 


is a simple example of a diagonal matrix. 


L0 0 lj 


Examples XIV 


l P = 

'i 2 r 

3 2 1 

Q = 

'3 2 f 

4 2 2 


.2 3 2. 


J 3 1. 


Find Q.P and P.Q. What does this illustrate? 
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1 

i 

2" 

1 = 

'1 

0 

O' 

2 

3 

2 


0 

1 

0 

J 

2 

1. 


_0 

0 

1. 


Find Bl. kl, B + kl, B-kl. 


"i 2 r 

Q = 

'1 0 1 

R = 

*i o -r 

1 0 1 


2 1 -2 


2 1 3 

.3 1 2. 


.3 2 3. 


J 2 3. 


Find PQ, ( PQ)R, QR. P(QR ). P(RQ\ K(P(?)- What conclusions 
can be drawn from these calculations? 


II 

' 3 

2 

r 

S = 

'2 

0 

O' 

T = 

‘2 

0 

O' 


-2 

3 

6 


0 

3 

0 


0 

-3 

0 


1 

-2 

4. 


_0 

0 

4. 


0 

0 

-4. 


Find AS, AT, SA, TA, S 2 , ST. 


II 

if) 

3 

2 

r 

V = 

'0 

0 

r 


-2 

3 

6 


0 

1 

0 


1 

-2 

4. 


J 

0 

0 . 


Find AV and VA. 


Geometrical Transformation of a (3 x 3) Matrix 

Just as (2x2) matrices cause transformations in geometrical 
figures of two dimensions, i.e. areas, so (3 x 3) matrices operating 
in three dimensions transform figures which are solids into new 
shapes. 






















































Introduction to Matrices 


Wc shall consider the unit cube using three mutually perpen¬ 
dicular axes ox, oy, oz, using a right-handed system as shown in 
the diagram. 

The points O. P, Q. R, S. T, U, V , defining the unit eube have 
the co-ordinates given in the table together with the co-ordinates 
of their images caused by the matrix transformation 


X 


'2 

■ nr*! 

or x' = 2x+y+z 

y 

— 

1 

2 T 

y' = x + 2y+2r 

_ z _ 


.1 

3 2.JLJ 

z' = x+3y+2z 


0 

p 

Q 

R 

S 

T 

U 

V 

0,0,0 

1,0,0 

0. 1,0 

0,0, 1 

1,0, 1 

0, 1, 1 

1,1,0 

1, 1.1 

O' 

F 

Q' 

R’ 

S' 

T 

U' 

V' 

0, 0,0 

2, 1, 1 

1,2,3 

1,2,2 

3, 3,3 

2,4,5 

3, 3,4 

4, 5,6 


These points when plotted in space give a figure whose faces 
are no longer squares but parallelograms and so the new solid 
is called a parallelepiped. The (3 x 3) matrix has transformed 
a unit cube into a parallelepiped of three dimensions. When 
considering the geometrical transformation caused by a (2x2) 
matrix we found that the area of the transformed figure and the 
area of the unit square were in the ratio equal to the value of the 
determinant of the matrix. When dealing with the (3 x 3) matrix 
the volume of the parallelepiped and the volume of the unit 
cube can be shown to be in the ratio equal to the determinant of the 
(3 x 3) matrix. 

If the transformed parallelepiped OP'Q'V'R'S'U'T is trans¬ 
formed by the second application of the matrix it is found that the 
new parallelepiped would have its volume increased again in the 
ratio of the value of the determinant since the solid OP'Q'V'R'S'U'T' 
could be considered as being composed of a number of unit volumes 
each of which would be increased in the ratio of the determinant 
and the new volume would be the sum of the parallelepipeds from 
each of these unit cubes. 
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The Calculation of the Determinant of a 
(3 x 3) Matrix 


Historically the idea of a determinant arose in the solution of 
linear equations and was found to be the value of the products 
resulting from the number of arrangements of the coefficients so 
that each coefficient came from a different row and column and was 
also given its appropriate positive or negative sign according to its 
position: we shall now proceed to demonstrate this. 

So with a (3 x 3) matrix the determinant formed from its elements 
must come from the arrangement of the elements taken three at a 
time since there are three rows and three columns. 


For the (3 x 3) matrix 


a x — b x —cv -Row 1 if wc select the 
a 2 b 2 c 2 


*>3 ^* 3 - 


Column 1 

element a x which is in column 1 and row 1 the elements not in its 
row or column are b 2 c 2 forming a (2 x 2) determinant. 

^3 C 3 

This (2x2) determinant is called the minor of element a x . To 
find the minor of element b x we strike out the elements in its row 1 
and column 2 leaving the elements a 2 c 2 1 to form its minor. 

C 3 


Every element of the matrix has its own minor and the deter¬ 
minant value of these minors can be calculated as shown earlier. 
Next if we give to each minor a positive or negative sign according 
to the position of the element of which it is the minor, shown by 
the pattern 

r+ - +■ 

- + - 

_+ - + 


then these signed-minors are called co-factors. 


In our matrix 

'a, b, c," 

the co-factor of a, is + 

b 2 c 2 


a 2 b 2 c 2 




no 

u 

no 

no 

_1 



= + (b 2 c$ b$c 2 

1 shown A x> 



73 




































Introduction to Matrices 


the co-factor of element />, is 
shown B u 

the co-factor of element c, is 
shown C,. 


a 2 c 2 
«3 C 3 


= -(fl 2 ‘'3-«3t'2) 




+ (a 2 b 2 —a i b 2 ) 


Definition of Determinant 

The value of the determinant is a,A + b l B + c 1 C which when 
multiplied out gives six products each formed from three elements: 


a,(b 2 c 3 -b 3 c 2 )+b,{ -( a 2 c 3 -a 3 c 2 )} +c l {a 2 b 3 - a 3 b 2 ) 


= a fiiC 3 + a 3 h x c 2 + a 2 b 3 c x -a x b 3 c 2 - a 2 b x c 3 - a 3 b 2 Ci 


Any row or column could have been used for the calculation 
of the determinant; had the first column been used, the value 
would have been 



b 2 c 2 

r 

6, c, 

) 

t> 1 c, 

a,A, +a 2 A 2 +a 3 A 3 = a, 


+ a 2 <- 



i>3 c 3 

L 

t>3 c 3 

J 

b 2 c 2 


= a^biCi-bjCi-a^b^i-biCj+a^biC^b^i) 

= a , b 2 c 3 +a 2 b 3 c l +a 3 b i c 2 -a x b 3 c 2 - a 2 b x c 3 - a 3 b 2 c x 

and this result is exactly equal to the previous value when the 
elements of the first row were chosen. The determinant when found 
is a pure number (a scalar). 


Example 

Find the determinant of the matrix 


1 

2 

3 


2 

4 

1 


1 

6 

2 . 


. The determinant 


is shown 


1 

2 

3 


2 

4 

1 


1 

6 

2 


using the same notation as with (2 x 2) matrices. 


(This notation was first introduced by Cayley in 1853.) 

Taking the elements of the first row 1 2 l their co-factors 


74 


are + 

4 6 

» 

2 6 

and + 

2 4 


1 2 


3 2 


3 1 


Section III 
or (8 — 6), —(4—18) and 


(2—12) i.e. +2, +14 and -10. 

Determinant (A) = (1 x +2)+2(+ I4) + (1 x —10) 

= 2 + 28-10 = +20 

The matrix used in the geometrical transformation on page 72 
shown in diagram 24 is f 2 1 1 

1 2 2 


2 1 1 
1 2 2 
1 3 2 


A = 2 

2 2 

-1 

1 2 

+ 1 

1 2 


3 2 


1 2 


1 3 


1 3 2. 

is the determinant of this matrix, call it A. 


note the negative sign of the 


second term. 

= 2(4-6)-1(2-2)+1(3-2) 

= -4+0+1 
= -3 

The volume of the parallelepiped is three times that of the unit 
cube (the negative sign has no significance with respect to volume 
but affects the position of the solid in space). 

When the determinant of a (3 x 3) matrix is zero the matrix is 
said to be singular and the transformation it represents causes 
the images of the eight points defining the unit cube to fall on a 
plane, i.e. the volume of the transformation is zero also and the 
images then define an area if only two row's or two columns are 


in the same ratio, as with the matrix 


The determinant 


is zero and rows 2 and 3 result in the y', z' of the transformed 
points being in a constant ratio of 1:2. All the points of this trans¬ 
form will lie in a plane which passes through the point O, O. O 
(the origin) and cuts the yz plane in a straight line which passes 
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through the origin and has the equation z = 2 y, this is shown in 
diagram 25. 

Since all the images of the points O, P, V, Q. R, S, U, T, lie in the 
plane shown, the figure has degenerated from a parallelepiped 
having volume, to a plane area. The singular matrix has resulted 
in an image of three-dimensional space in two-dimensional 
space. This topic will be dealt with in another aspect later, when 
the subject of Linear Dependence is studied. 



Although we have shown that the operation represented by a 
(3 x 3) matrix has a geometrical meaning wc can consider (3 x 3) 
matrices as an independent system having an algebra of their own 
as was done earlier with (2x2) matrices. Wc now proceed to a 
study of these matrices independent of the operations which they 
represent. 


Examples XV 

1 Find the value of the determinant |P| = 


1 

3 

2 


2 1 
2 1 

3 1 


of the 
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matrix P — 


1 

3 

2 


2 I 
2 1 

3 2. 


2 Find the value of the determinant |(?| of the matrix 

q = r 3 2 r 

4 2 2 
1 3 1 


3 Find the value of the determinant |/l| of the matrix 


A = 


1 0 3 


2 1 4 
.0 1 3 


4 


Is the matrix B = 


'2 1 1 
1 2 2 
_2 4 4 


singular or non-singular? (Hint: 


find the value of its determinant |B|.) 


5 Using the matrices A and B of questions 3 and 4, find C = A. B, 
then find the determinant of matrix C. Is matrix C singular or 
non-singular? 


6 Find the determinant of 


7 Find the determinant of 


2 0 O' 
0 2 0 
.0 0 2 . 
2 0 0 ‘ 
0 3 0 
0 0 4 


8 Find the determinant of D = 

'1 0 O' 

and £ = 

'1 0 O' 


2 2 0 


0 2 0 


4 3 2. 


O 

o 

to 
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The Transpose of a Matrix 

The transpose of a matrix A is that matrix whose rows are 
the columns of A and whose columns are the rows of A and is 
denoted by A or A 



*1 

c,- 


~ Cl\ 

a 2 


N 





\ 

«2 


<■2 

A' = 

b x 

\ 



\ 




\ 


-«3 

b 3 

\ 

c 


Sx 

c 2 



Inspection shows that the elements along the leading 
diagonal are unchanged in position and since we can calculate 
the value of the determinant using the elements of a row or the 
elements of a column then the determinant value must remain 
unaltered. 

If we denote the determinant of the matrix A by the symbol 
\A\ then 

\A\ = \A‘\ 

If the transpose of a matrix gives a matrix which is equal to the 
original, such a matrix is said to be symmetric. 

Example 


Q = 

'1 2 3' 

Q' = 

“ 1 2 3“ 


2 4 1 


2 4 1 


.3 1 2. 


.3 1 2. 


Q = Q' 

Q is symmetric. 


The Adjoint Matrix 


Definition 

The matrix formed by transposing the elements formed from 
the co-factors is called the adjoint matrix. 


For the matrix A = 


3 1 2 


used in the example in page 74 
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the co-factors of the elements are by calculation of the minors 
2 14 -10' 

-3-1 5 

8-4 0. 

The Adjoint matrix is formed by transposing the rows and 
columns. 

Adjoint matrix is 2—3 8’ 

14 -1 -4 

.-10 5 0J 

Now this Adjoint matrix has very special properties in relation 
to the matrix A from which it was formed and the Unit matrix /. 
This will be shown by multiplication. 


A .[Adjoint A\ = 


1 2 1 
2 4 6 
.3 1 2. 

20 0 
0 20 0 
. 0 0 20 J 


' 2 

14 

.-10 
0^ 


-3 

-1 

5 


8' 

-4 

0 


= 20 


1 0 0 


A. 


f Adjoint A 


= 20 / 


1- 


L" 20 

But on page 75 we showed the determinant of A to be 20 

[ Adjoint A 

- 


(0 


I/ll 

Next we find 


= / 


iriH-i 


1 

2 

-3 

8' 

“1 

2 

r 

20 

14 

-1 

-4 

2 

4 

6 


_—10 

5 

0. 

.3 

1 

2. 
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1 

20 


20 
0 
L 
1 0 
0 1 
L0 0 


= / 


0 

20 

0 

O' 

0 

1. 


0 

0 

20 . 


... (ii) 


n r- ~ \ = 1 ( sh ° wn at (0) 

and |^J = / ( s h own at («)) 

The matrix [ Adjoint A~\ can also be written 1 [Adjoint A 

L Ml J Mir 


since the determinant \A\ is a number, i.e. a scalar, and this matrix 
has special properties since pre-multiplication or post-multiplication 
of A by this matrix results in the unit matrix. It is unique and is 
the inverse of A, i.e. A~ l . 


A~ l 


Adjoint A ] 


If for any (3 x 3) matrix the determinant is zero then the inverse 
of that matrix cannot be defined. Such a matrix has already been 
termed a singular matrix. But for all non-singular matrices the 
inverse matrix can always be calculated from the Adjoint matrix 
as shown. [The method holds for all non-singular square matrices 
of higher orders also—it is a general method for square matrices 
of any order but at this stage we are only concerned with (3 x 3) 
matrices.] 


Examples XVI 

1 Find the transpose A'of A 


'1 2 3' 
3 2 1 
.2 3 2. 
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2 Find the transpose B' (or B ) of B 


3 2 1 

-2 3 6 

1 -2 4. 


3 In questions 1 and 2 compare the leading diagonal of A' with 
that of A and of B' with that of B. What conclusion do you draw? 


4 Form the matrix C whose 
the elements of Q = T3 2 1 

4 2 2 

1 3 1. 


elements are the co-factors of 
, then by transposing the elements 


of C find [adjoint Q]. Calculate the determinant of Q, i.e. |Q|. 
Hence state Q ' 1 and verify that QQ~ l = Q' l Q = I- 

5 Find A ~ 1 and verify A~‘.A = A.A~' = / when 


A = 


6 If P = 


7 If R = 


8 If A = 


'a h g 
h b f 
.g f c_ . 
notice about A' and A? 


1 0 3' 

2 1 4 
L0 1 3. 

1 2 1 " 

3 2 1 
2 3 2. 

-1 2 O' 

-3 1 2 

3 0 -4. 

write out the transpose of A. What do you 


find P 1 and verify. 


find R 1 and verify. 


9 If P = 

" 1 

2 

f 

Q = 

'l 

1 

r 


1 

0 

1 


2 

1 

2 


.3 

1 

2. 


.3 

2 

3 


find PQ, P\ Q\ QP , 
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P'Q\ Q'P\ ( PQY , (QP)' and show that 


(PQY = Q'P 
(QP)' = P'Q'. 

(P‘ = P' Alternative symbols). 


10 If P = 

"1 

2 

1 

, Q = 

l 

1 

r 

, R = 

'1 

0 

r 


1 

0 

1 


2 

1 

2 


2 

1 

3 


.3 

1 

2. 


.3 

2 

3. 


.1 

2 

3. 


PQR, QR, P\ (QR)'P! R’Q', R'Q'P' and hence show that 

( PQR)' = (QR)'P' 

= R'Q'P’ ... (the Reversal Rule) 
(P' = P' P' means P') 


The Solution of Three Simultaneous Equations 
Using Matrix Method 


if 


—x + 2y—3z = -8 
2x—y + 4z = 17 
3x + 4y +z = 22 


then 


‘-I 

2 

-3“ 

X 


"-8' 

2 

-1 

4 

y 

= 

17 

- 3 

4 

1. 

_Z 


. 22. 


... (0 


X 


"-8" 

y 

= 

17 

7 


.22. 


If the 3 x 3 matrix is represented by A. 


We now calculate the inverse of A in three steps, 
(i) Determinant of 


A = \A\ = — 1( — 1 — 16) —2(2— 12)—3(8 + 3) = 4. 
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(ii) Adj.A = 


'-17 

10 

11" 



-14 

8 

10 



. 5 

-2 

-3. 

(iii) A~ 1 = 

‘ Adj. A 

™ J 

1 

" 4 

'-17 





-10 
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'-17 

-14 

5 

10 

8 

_ 2 

. 11 

10 

-3 


-14 

5 ' 

8 

_2 

10 

-3. 


We now pre-multiply both sides of the matrix equation (i) 


-17 -14 


10 

11 


8 -2 
10 -3 


-1 2 -3' 

2-1 4 

. 341 . 
'4 0 0 
0 4 0 
.0 0 4. 
1 0 0 
0 1 0 
.0 0 1 . 


-17 -14 


10 

11 


8 

12 

16 


8 -2 



— 8 ' 


17 


_ 22 . 


For these two column matrices to be equal, their corresponding 
elements must be equal 

x = 2 
=> y = 3 
2 = 4) 


The result can be checked by substitution or by working in the 
normal algebraic method. We would not normally use the matrix 
method for this kind of problem, but at a later stage where machine 
or computer facilities are available similar applications are possible. 
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Introduction to Matrices 

Manufacturers of desk calculating machines often supply in¬ 
structions with the machines for this type of matrix solution and 
hold short courses at many centres on the subject. With large 
electronic computers other methods are available when equations 
with hundreds of unknowns have to be solved, but these methods 
are only understood after matrix methods have been mastered. 


Examples XVII 


in matrix form P 

x~ 

= 

0" 


y 


2 


_z_ 


.2. 


x + 2y + z = 0 
3x + 2y + z = 2 
2x + 3y + 2z = 2 

then find P~ l and solve the equations using the matrix method. 
2 Express the equations 

-x + 2y = 5 
— 3x + y + 2z = -1 
3x -4z = +9 

in matrix form and find the solutions using matrix methods. 


3 If Q = 


by matrix methods. 
4 Solve 


"3 2 r 

and Q 

x" 

= 

16' 

4 2 2 


y 


22 

.1 3 1_ 


z_ 


.15. 


, find the set of solutions 


x + 3z = 7 ^ 
2x + y + 4z = 13 
2y+3z = 12 j 


by matrix methods. 
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Product of Two Matrices and their 
Determinants 


On page 74 matrix A = 


"I 2 1" 

2 4 6 
.3 1 2. 


has a determinant value 


of 


20 (found on page 75), we call the matrix B = 


-1 

2 


which was used on page 82 and has determinant 4, 
and \B\ = 4. 

Mpvi \x tt* find A R z inH R A ' 


2 -3' 

-1 4 

4 1_ 

i.e. |/4| = 20 



'1 

2 

1 

-1 

2 

-3' 


' 6 

4 

6' 

A.B = 

2 

4 

6 

2 

-1 

4 

= 

24 

24 

16 


.3 

1 

2. 

L 3 

4 

1 . 


. 5 

13 

-3. 


If we call this matrix C such that A.B = C we now find that 
|C| = 80 but M| = 20 and |B| = 4. 

=* MI-MI = |C| although matrix C is a new matrix obtained by 
pre-multiplying B by A. 

Next we find the product B.A = D 



"-1 

2 

-3' 

‘1 

2 

f 


'-6 

3 

5" 

B.A = 

2 

-1 

4 

2 

4 

6 

= 

12 

4 

4 


. 3 

4 

1 . 

.3 

1 

2. 


. 14 

23 

29. 


The matrix D is obtained by post-multiplying B by A and is 
obviously different from C. 

Once more we have shown that A.B # B.A. 

But although multiplication of matrices is non-commutative 
since the value of |D| = 80, |B| MI = Ml- Generally, 

\A\.\B\ = \A.B\ 

and 


|B|.MI = \B-A\ = \A.B\. 
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Thus the multiplication of the determinants does follow the 
commutative law and this must be so, because \A\, |B|, |D|, |C| are 
all scalar numbers and so follow the laws of number algebra. 

In the algebra of (2 x 2) matrices it was shown that the product 
of any matrix with a singular matrix gave a singular matrix 
as the product. We have now demonstrated that the product of 
two (3 x 3) matrices has a determinant which is the product of the 
determinants of the two separate matrices. This is always true 
and hence with all matrices the product of a matrix with a singular 
matrix is always a singular matrix. This result is important and 
has many implications in later work and we shall use this property 
when dealing with the characteristic equation of a matrix, which 
has many applications especially in quantum mechanics. 


The Calculation of Determinants of the Third 
Order 


It was shown earlier that for the third order determinant 


a i b i 
a 2 b 2 

a 3 b 3 


Cl 

C 2 

C 3 


the value of the determinant A is a l A l + b l B l +CiC, where 
A t , B,, and C, are the co-factors of a u b t and c,. 


A = a l A i +b l B l +c,C, 


= a, 

b 2 c 2 

-b t 

a 2 c 2 

+<T 

a 2 b 2 


b 3 c 3 


c 3 


a 3 b 3 


= (•i(b 2 c i -b 3 c 2 )-b l (a 2 c 3 -a 3 c 2 )+c l (a 2 b 3 -a 3 h 2 ) 

= Oib 2 c 3 — a x b 3 c 2 — a 2 b i c 3 + u 3 bit' 2 +a 2 b 3 c | —a 3 b 2 c l 


= a l b 2 c 3 + a 2 b 3 c t + a 3 bic 2 -a l b 3 c 2 -a 2 b l c 3 -a 3 b 2 c l 


There are three positive products and three negative products 
and it can be seen that they all have their factors lying along 
diagonals drawn through the array of elements. 
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x / x 

a 3 b 2 C i -a } b 3 c 2 -U 2 frl c 3 



+ a i b 2 c 3 +a 3 b,C 2 


+ n 2 b 3 Ci 


This provides a rapid way to evaluate the positive and negative 
products—usually it can be done mentally. It is known as the 
rule of sarrus but is only applicable to third order deter¬ 
minants. 


Examples XVIII 

I Use the Rule of Sarrus to evaluate 


and 300 
0 4 0 
0 0 2 

Find the determinants of the matrices A = 


1 

0 

3 

and 

i 

0 

2 

2 

1 

4 


3 

4 

5 

0 

1 

3 


5 

6 

7 


TOO 
2 3 5 


and 


B = DO 6 
3 4 15 
.5 3 20. 


3 Is matrix M = 


4 1 3J 

. Then multiply A and B and show that 


\A.B\ = \A\.\B\. 


0 2 3 

-2 0 4 

L—3 -4 0J 


singular or non-singular? 


4 If matrix N = 

'2 -2 2' 

and M = 

' 0 2 3' 


1 3 2 


-2 0 4 


. 2 44 3. 


O 

1 

1 

_1 


is the 


product N .M singular or non-singular? 
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5 Use the Rule of Sarrus to show that 


U’l+di) 


n, fc, c, 


a, h, d, 

a 2 b 2 (^* 2 +^ 2 ) 

= 

"2 b 2 c 2 

+ 

a 2 b 2 d 2 

a s ^3 (0 + ^ 3 ) 


b 3 c 3 


a 3 b 2 d 3 


II 

NC 

1 

0 

0 

, find det T. i.e. |T|, and U = 

1 

5 

7 ], 


3 

1 

0 


0 

1 

3 


4 

2 

1. 


.0 

0 

L 


find det U , i.e. \U\. 

1 If A = TU where T and V have the values of question 6 
without calculation give \A\. 


Reflections in Three Dimensions 


The diagonal matrix P = 


1 0 0 

0 1 0 

0 0 - 1 . 


When treating 2x2 matrices it was shown that the matrix 
1 0 represents a reflection about the x axis. The matrix P 

0 - 1 . 

above, is similar in form and if we examine the transformation 


given by 


1 

0 

0' 

x" 

= 

X 

0 

1 

0 

y 


y' 

0 

0 

-L 

_Z_ 


.z'J 


Multiplication gives 



The image T has the same ordinates in the xy plane but is 
obviously a reflection in the xy plane. 
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The matrix 


plane. 

The matrix 


T 

0 

0 


0 

-1 

0 


0 

0 

1 


represents a reflection in the (x, z) 



0 O' 
-1 0 
0 1 . 


represents a reflection in the (y, z) 


plane followed by another reflection in the (x, z) plane. 

It is obvious that all the possible reflections are included in the 
eight (3 x 3) matrices 



'1 

O 

O 

'l 

0 

O' 

I 

0 

O' 


-1 

0 

0 ' 



0 

1 0 

0 

1 

0 

0 - 

1 

0 


0 

1 

0 



.0 

0 

1. 

.0 

0 - 

- 1 . 

.0 

0 

1 . 


0 

0 

1 . 


1 

0 

O' 

~ 

-1 

0 

O' 

'-1 

0 

0 " 


1 

0 

O' 

0 - 

-1 

0 


0 

-1 

0 

0 

1 

0 


0 - 

-1 

0 

,0 

0 

-1. 


0 

0 

1 . 

L 0 

0 

— 

1 . 


0 

0 

- 1 . 


The eight possible reflection matrices form a group of possible 
reflections in the three planes (x, y), (y, z), (x, z). 


D 
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The matrix 


'2 0 0 ' 

0 1 0 

.0 0 - 1 . 


therefore represents a reflection in 


the (x, y ) plane and a stretch of the x ordinate. 


The matrix 


would represent a reflection in the 


LO 0 -1J 

(x, y) plane and a stretch in two directions in that plane. 


The Diagonal Matrix 


Some of the matrices studied so far, the Unit matrix 


1 

0 

0 


o 01 
1 0 
0 1. 


the Scalar matrix 


k 0 0 
0 k 0 
0 0k 


= kl, and the Reflection matrices 


are sub-groups of a very large and important group of matrices 


known as diagonal matrices, of the form 


a 0 0 


where 


0 p 0 
0 0 y. 


a, fi and y are scalars. Using the method of expansion by co¬ 
factors or by using the Rule of Sarrus, it is obvious that the 
determinant of the diagonal matrix above is a/?y. Since this is the 
general diagonal matrix, we can state in general terms that the 
determinant of a diagonal matrix is the scalar obtained by multi¬ 
plying together the elements along the diagonal. 

Diagonal matrices play a very important part in later work and 
require some special treatment. 

Now take two different diagonal matrices 




0 

0 ' 



"“2 

0 

0 ■ 

A = 

0 

bx 

0 

and 

B = 

0 

b 2 

0 


.0 

0 

Cl. 



.0 

0 

Cl. 
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0 

O' 

a 2 

0 

0 " 


'a,a 2 

0 

0 ' 

A.B = 

0 

bx 

0 

0 

b 2 

0 

= 

0 

bxb 2 


also 

.0 

0 

C|_ 

.0 

0 

C 2- 


0 

0 

CxCl. 


‘“2 

0 

0' 

'a, 

0 

0' 


a i a \ 

0 

0 ' 

B.A = 

0 

bi 

0 

0 

bx 

0 

= 

0 

b 2 b, 

0 


.0 

0 

c 2- 

0 

0 

Cx- 


L o 

0 

C 2 C 1- 


Hence these two diagonal matrices commute 


A.B = B.A 


and the product is also a diagonal matrix. We induce that 
A. B. C - C. B.A if C is a diagonal matrix the product must be 
another diagonal matrix. 

It follows too that if A = B then 


and hence 



'a? 

0 

1 

O 


'ai 

0 

0 

A 2 = 

0 

bi 

0 

II 

0 

bi 

0 


.0 

0 

r 2 
Ci J 


.0 

0 

Cl- 


al 


A" 


0 

.0 


0 0 ' 

IS[ 0 
0 d[. 


Diagonal matrices are very simple to use and because of their 
simplicity are useful in many applications. The inverse is simple 
to And. 





~bxt 

■1 0 

0 ' 

A 

-i _ 

1 

0 

o,c, 

0 



0,6,0, 







. 0 

0 

axbx. 


a~ 1 

0 

0 ' 



= 

0 

6 -* 

0 




. 0 

0 < 

.-1 




This latter result is easy to remember. 
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Triangular (3 x 3) Matrices 

Triangular (3x3) matrices have properties similar to those of 


(2 x 2) matrices. If A = 

'a, 

by 

<T 

and T = 

T 

1 

O' 


a 2 

b 2 

<2 


0 

1 

0 


-«3 

b 3 

<3- 


.0 

0 

1. 


triangular matrix because the elements below the leading diagonal 
are all zeros. 



'1 

1 

0 


'ay 

b i 

*r 


'a, + a 2 

b 1 +b 2 

Cy+c 2 

T.A = 

0 

1 

0 


«2 

b 2 

C 2 

= 

a 2 

b 2 

C 2 


.0 

0 

1 . 


J>3 

b 3 

c 3 _ 


- 

b 3 

C 3 - 


Pre-multiplication of A by T has been equivalent to the addition 
of row 2 to row 1. 



'01 

by 

Cy~ 

1 

1 

O' 


’<*1 

a, +6, 

Cl' 

A.T = 

02 

b 2 

f 2 

0 

1 

0 

= 


a 2 -f* b 2 

Cl 


.03 

b 3 

c 3 - 

.0 

0 

L 


-a 3 

03 + ^3 

C 3- 


Post-multiplication of A by T has caused the addition of column 1 
to column 2. 



'1 

0 

3' 



T 

0 

3' 

'fli 


Cy 

If T, = 

0 

1 

0 

then 

TyA = 

0 

1 

0 

«2 


Cl 


.0 

0 

I. 



.0 

0 

1. 

-«3 

b 3 

c 3 . 


0,4-303 b!+36 3 c,+3c 3 ’ 

= o 2 b 2 c 2 
- 03 *>3 c 3 . 

The matrix T x has resulted in the addition of three times the 
third row to the first row. 

We can state that pre-multiplication by a suitable triangular 
matrix causes the addition (or subtraction) of one row (or a 
scalar multiple of that row) to another row. Post-multiplication 
has the same effect on columns. 
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If N = 

‘0 

1 

O' 

and as before A = 

'a , 

by 

Cy 


1 

0 

0 


«2 

b 2 

C 2 


_0 

0 

1. 


-a 3 

b 3 

C 3 . 



'0 

1 

O' 

'ay 

by 

Cy' 


~a 2 

b 2 

c{ 

N.A = 

1 

0 

0 

a 2 

b 2 

c 2 

= 

ay 

by 

C 1 


.0 

0 

I. 

a 3 

b 3 

C 3- 


.a 3 

b 3 

^3- 


Pre-multiplication by N has interchanged row 2 and row 1 of 
matrix A. (Post-multiplication would have interchanged column 1 
and column 2.) 


The diagonal matrix M = 


1 0 0 " 
0 3 0 
0 0 2 . 


would give a result as 


follows: 



'1 

0 

O' 

’*1 

by 

Cy' 


'a y 

by 

Cy ' 

M.A = 

0 

3 

0 

a 2 

bi 

c 2 

= 

3q 2 

3 b 2 

3 c 2 


.0 

0 

2. 

a 3 

b 3 

c 3 . 


.2 a 3 

2b 3 

2f 3- 


The second row has been multiplied by 3 and the third row by 2. 


Examples XIX 


1 A = 

'3 -1 2' 

1 1 0 

. D = 

'2 0 0" 

0 3 0 

, find AD and DA. 


.2 -2 0. 


.0 0 -2. 


2 B — 

T 2 31, 

D = ' 

1 0 0T find BD and DB. 


2 0 4 
12 6 


0 4 0 
L0 0 7. 


3 E = 

'I 

0 

0" 

. F = 

'3 

0 

0 


0 

-1 

0 


0 

2 

0 


.0 

0 

3. 


.0 

0 

-1. 


, find |£| # |F|, \EF\. 
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Introduction to Matrices 


4 G = T2 0 0 

0-2 0 

L0 0 3. 

5 B = T3 — 1 2], T = 

1 1 2 

.3 -2 0 

then find |B|. |7|, |BT|, |TB|. 


, find |G|, |G J |, |3G|. 


110 
0 1 0 
0 0 1. 


find BT and TB. 


6 

T = 

'1 

1 

O’ 

, y = 

'1 

0 

O’ 



0 

1 

0 


1 

1 

0 



.0 

0 

1 . 


.0 

0 

1 . 

7 

L = 

1 

0 

O’ 

, u = 

’1 

d 

e~ 



a 

1 

0 


0 

1 

f 



h 

c 

1 . 


.0 

0 

1 . 


, find TV, VT, V 2 , T 2 


, find the matrix A such that 


A = LU. To which class of matrices do L and V belong? 


8 If in question 1,A = 


'1 2 2 
2 5 8 


, find the values of a, b, c, d, e,f 


.3 7 11. 

in L and U. Express A as the product of two triangular 
matrices. 

9 Find the inverse of matrix L in question 8. 

10 Find |L| and |l/| for the matrices of question 8. 


The Elementary Transformations 

We have shown that three elementary transformations 
of a matrix can be brought about by pre-multiplication or post¬ 
multiplication by another suitable matrix. 

Definition 

The elementary transformations of a matrix are 
(i) the interchange of two row's (or columns) 
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(ii) the multiplication of each element of a row (or column) by 
scalar 

(iii) the addition of a scalar multiple of one row (or column) to 
another row. 

If matrix A can be transformed by one of these elementary 
transformations into matrix B, then A is said to be equivalent 
to B. (Shown A — B.) The elementary transformations can be 
used to find the inverse of a matrix. 


Finding the Inverse Matrix using the 
Elementary Transformations 

The work is carried out in two parallel columns and since each 
transformation is the result of multiplying by a suitable matrix, 
each matrix operation can be represented by a different letter. 
Pre-multiplication causing row operations will be used on the unit 
matrix also. 



*i 2 r 



O 

o 

A = 

2 4 6 


/ = 

0 1 0 


.3 1 2. 



o 

© l 


Subtract twice row 1 from row 2. This is equivalent to pre¬ 
multiplication by some matrix P. 


Both matrix A and matrix / are subjected to the same matrix 
operation on each occasion. 


PA = 

'12 1' 

0 0 4 

.3 1 2. 


PI = 

'10 0' 

-2 1 0 

.001. 


Subtract 2 

QPA = 

1 x row I from row 3 

'i 2 r 

0 0 4 

.0 -5 -1. 

QPI = 

'10 0' 

-2 1 0 

.-3 0 1. 
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Introduction to Matrices 
Add row 3 to row 1 



1 -3 O' 



'-2 0 1* 

RQPA = 

0 0 4 


RQPI = 

-2 1 0 


.0 -5 -1. 



.-3 0 1. 


Add 4 x row 3 to row'i 



'1 -3 O' 



" -2 0 r 

SRQPA= 

0 -20 0 


SRQPI = 

-14 1 4 


_0 -5 -1. 



.-3 0 1. 


Subtract | x row 2 from row 3 



"1 -3 O' 



" -2 or 

TSRQPA = 

0 -20 0 

.0 0-1. 


TSRQPI = 

1 - 

1 

M*- 4^. 

1 

— 

O ^ 

I. 


Divide row 2 by — 20 and 
row 3 by — 1 



1 

-3 

0 



'-2 

0 

1 ' 

UTSRQPA = 

0 

1 

0 


UTSRQPI = 

7 

10 

-1 

20 

-1 

5 


.0 

0 

1. 



1 

1 

0 






~2 

4 


Add 3 x row 2 to row 1 


WVTSRQPA = 


1 

0 


0 


0 

1 

0 


0 


0 

I 


WUTSRQPl = 


10 

1_ 

10 

~2 


-3 2 
20 5 

-1 J 
20 '5 
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In column 1 the matrix A has been transformed to the unit 
matrix 


(WUTSRQP).A = / 

** (WUTSRQP) = A' 1 

But in the set jnd column we have (WUTSRQP).I = A~'I = A *. 
Thus the matrix in column 2 is the inverse matrix A -1 . 


• 

1 

10 

-3 

20 

2' 

5 


2 

-3 

8" 

A~ l = 

7 

10 

-1 

20 

-1 

5 

1 

20 

14 

-1 

-4 


-1 
. 2 

1 

4 

0 


-10 

5 

0. 


The process has been set out in full, since the operations are 
simple many of them can be carried out in a single operation 
and the result then follows quickly. It must be emphasized that 
row operations only have been used, i.e. pre-multiplication. 

Although we have illustrated the elementary transformations 
using (3 x 3) matrices, they can be extended to other groups of 
matrices. 

Further applications of the elementary transformations arise 
in later work. 

For finding the inverse of matrices larger than (3 x 3) this method 
is more efficient than any other method using determinants. 


Examples XX 

1 Using elementary transformation by rows reduce 
the inverse matrix and check. 


3 4 
2 3 


to 


2 Find the determinant of 


-2 


and then reduce to the 


4 -7 

inverse matrix using row transformations. 

3 Using elementary row transformations find the inverse of 


2 1 1 
1 2 2 
1 3 2 


. Check your result. 
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Introduction to Matrices 


4 A = 


5 Q = 


6 R = 


1 0 
2 1 
0 1 
3 2 


, find A 


- i 


, find Q 


- i 


0 1 


0 1 0 

Li o o. 


, find R 


-l 


Applications of Matrix Algebra 

Matrix Algebra is one of the great developments of nineteenth- 
century mathematics but it is only in recent times that its applica¬ 
tion in so many different fields has been realized. The rapid 
growth of Linear Programming since 1947 has emphasized the 
need for matrix methods in this new branch of mathematics. 
More recently the application of matrix algebra has given the 
electrical engineer a powerful method in the analysis of electrical 
networks. 

The use of matrix methods in electronics, quantum theory, 
newer branches of social science, aerodynamics, theory of vib¬ 
rations, physical chemistry, higher mathematics and above all 
in modern computer techniques, makes certain concepts of matrix 
algebra essential for the understanding of these applications. 
Of outstanding importance for this purpose is the study of Eigen¬ 
values and Characteristic Function, Linear Dependence and Rank. 
Because of their wide applications these and other topics will be 
dealt with at appropriate places in the succeeding chapters. 


Latent Roots, Eigenvalues and Characteristic 
Equation 


The (2 x 2) matrix 



4 

13 


represents a transformation in 


two-dimensional space. 
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So a vector 

x' 

is transformed to 

x' 

, i.e. A. 

X 

= 

x' 






.y . 




If the transformed vector 


y J 


has the same direction as the original 


vector 


then we have shown earlier that it can be represented by 


iyl 


kx 

or a 

x‘ 

.'-y . 


.y. 


Such a vector is called a characteristic 


vector (or, eigenvector) of the matrix A. 


In the case of the (2x2) matrix there are not more than two 
values of /. possible. The values of k are called the charac¬ 
teristic ROOTS or LATENT ROOTS Or EIGENVALUES of the 
matrix. We shall also show that there are not more than three 
characteristic roots for a (3 x 3) matrix and so on. 

If we use X to denote the vector then A. X = kX for the condition 
understood, and since kX = klX where / is the Unit Matrix then 


AX = klX 

(A-kl).X = 0 ... (i) 

The determinant of the R.H.S. is zero, therefore one of the two 
matrices (A-kl) or X must be singular and since X is a vector 
then (A-kl) must be singular, i.e. its determinant is zero. 


\A-kI\ = 0 


This is called the char acteristic equation of the matrix A. 
\(A-kl)\ is called the characteristic function of A. 

The solution of \A —kl | = 0 gives the latent roots or eigenvalues. 
Substituting the known values of the elements of A gives 


[2 4' 


ol 

1.3 13. 


.o ;J 


2-k 4 

3 13-;. 


= 0 

= 0 
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Introduction to Matrices 


- (2 —;.)(13—A)—12 = 0 
A 2 — 15A + 14 = 0 
- (A-14XA-1) = 0 


=* A = 14 or l 

Substituting these values for A in the equation shown at (i) gives 

X = 0 


[2 

4 


14 

oj 

b 

13. 


. 0 

14 J 


-12 4 

3 -1 


X = 0 


If vector X is 


then 


-12 4' 

X 

. 3 -1. 

.y. 


= 0 



-* — 1 2x + 4y = O') 

x = 1| is an c 

and 

3 x-y = oj 

y = 3j 


For A = 

1 the relation becomes 

'2-1 4 

x" 



. 3 13-1. 

.y. 


>x+4)' = 0 


and 


3x + 1 2y = 0 


y 

= -1 

The two vectors 

T 

.3. 

and 

4' 

. -1 . 


= 0 

x = 4] satisfies these equations 

are the eigenvectors of the 


matrix A. 


Any vector of the form k 

r 

or k' 

' 4' 


.3. 


. - 1 . 


, where k and k' are any 


real numbers, would also satisfy the relation (A-/.I)X = 0. 
So finally we state that any vectors whose elements are in the ratio 
of 1:3 or 4: — 1 are eigenvectors of the matrix A. 


Section III 


Next we apply the transformation A to them 


'2 

4' 

T 


X 1 ' 

i 

II 

£ 

II 

14 

T 

.3 

13. 

.3. 



<N 

II 


.3. 


'2 

4' 

4' 


' x"" 

II 

II 

** 

It 

4' 

.3 

13. 

.-1. 



y" = -1 

.- 1 . 



The vector OP has transformed to OF in the same direction but 
14 times magnified in the first case and magnified by 1 in the second 
case when OQ transforms to OQ in the same direction. 


Any vectors of the form k 


or k! 


will be transformed by 


the matrix A in the same direction. No vectors of 2 dimensions in 
any other direction will be transformed by matrix A in the same 


direction. 

We have given a geometrical 
of the matrix A = 

3 13 


illustration of the eigenvectors 
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Introduction to Matrices 


Reduction of a Matrix to Diagonal Form 


The matrix constructed from the two eigenvectors 

T 

and 

4' 


.3. 


-1. 


has important and useful applications. Let this matrix be 
Q = T 1 41, its inverse 


1 4 
3 -1 


Q~ l = 


-1 

13 




It can be proved that if the original matrix A is pre-multiplied by 
Q 1 and post-multiplied by Q the new matrix will be a diagonal 
matrix, call it N 


N = Q~ l AQ = 


-1 

-1 

-4' 

'2 

4' 

1 

4' 

13 

.-3 

1 . 

.3 

13. 

.3 

- 1 . 



f 14 

4' 

. [42 

- 1 . 


_ i — 182 O' 
13 . 0 -13 


14 0 
. 0 1 . 

Note that the elements along the diagonal are the original roots 
of the characteristic equation of matrix A. 


Examples XXI 

1 Write down the characteristic equation and solve for the 

latent roots (or eigenroots or eigenvalues) of ,4 = [-15 -12 

. 24 19_ 

2 Write down the characteristic equation and find the latent 

roots (or eigenroots) of B = [ 3 — 1' 

0 4 
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3 Find the latent roots of C 


17 24 

-12 -17 


4 Find the eigenvectors of matrix A in question 1. 

5 Form the matrix Q with the eigenvectors of question 4, find 
its inverse and hence find the diagonal matrix equivalent to matrix 
A of question 1. 


6 


Reduce matrix B = 


3 

0 


-1 


to diagonal form. 


4 


Latent Roots and Eigen Vectors of a (3 x 3) 
Matrix 


To find the latent roots and eigenvectors 


2 3 

1 1 


— 2 " 
2 
3. 


of the (3 x 3) matrix 


The characteristic equation is |J5 — A/| = 0, from which the 
latent roots can be calculated: 



O 

1 

1 


1 

O 

O 



2 3 2 

- 

O 

O 



|- 

U> 


1 

O 

O 



(o-;.) -l 


2 

1 


(3-A) 

1 


-2 

2 


= 0 


(3-A) 


A 3 — 6/ 2 + 1 1a — 6 — 0 
^ (A-1XA—2XA-3) = 0 
=»/.= ! or 2 or 3 

The latent roots of matrix B are -I-1, +2, and +3. 
Substituting these roots in turn in the matrix relations 




V 


~x~ 


V 

y 

= 0, [B—(2)/] 

y 

= 0, [B—(3)/] 


z_ 






= 0 
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Introduction to Matrices 
gives 

(i) 


'-1 

-1 

-2" 

V 

2 

2 

2 


1 

1 

2 

_z_ 


= 0 


—x—y—2z = 0 
2x + 2 y+z = 0 
x+y+2z = 0 
x = 1 , y = — 1 , z = 0 


The eigenvector corresponding to the root + 1 is 


r 

-1 

0 


(ii) 


'-2 

-1 

-2" 

V 

2 

1 

2 

y 

. 1 

1 

1, 

_z_ 


= 0 (for root X = +2) 


-* 2x+y+2z = 0 
2x+y + 2z = 0 
x+y+z = 0 

y = 0, x = + 1 , z = - 1 
The eigenvector corresponding to root + 2 is 


r 

0 

- 1 . 


(iii) 


'-3 

-1 

—2 

x" 

2 

0 

2 

y 

. 1 

1 

0. 

_z_ 


= 0 (for root / = +3) 


3x+.y+2z = 0* 

2x + 2z = 0 
x + y = 0 

x = 1, y = -1, z = - 1 
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The eigenvector corresponding to root + 3 is 1 

1 

-1 

The matrix Q formed from these three eigenvectors is 

" i i r 

e= -1 o-i 

0-1-1. 

and if the inverse of Q is calculated, it is found to be 

'-1 0—11 

Q-' = -1 


-1 -1 
1 1 


Then 


Q~ l BQ = -1 


-1 0 - 1 ' 
-1 -1 0 


1 1 ULl 


0 -1 -2 
2 3 2 


1 



" 1 

l 

r 


-1 

0 

-l 


0 

-i 

-l. 




= -1 


-1 0 -f 

-1 -1 0 
111 . 

-1 0 0 
0-2 0 
0 0 -3J 


1 2 3 

-1 0 -3 

0 -2 — 3J 


1 0 0 

0 2 0 

0 0 3. 


= M 


The operation has produced a diagonal matrix with the three 
eigenroots (or latent roots) of B as elements along the leading 
diagonal. 
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Introduction to Matrices 


Examples XXII 


1 


If A = 


— 2 -3 -1 


write down the characteristic equation 


1 2 1 

. 3 3 2. 

and solve for the latent roots. 


2 Find the eigenvectors for the matrix A above and hence reduce 
A to diagonal form. 


3 If B = 


1 -1 -2 
2 4 2 


write down the characteristic equation 


1 1 4 J 

and solve for the latent roots. 


4 Find the eigenvectors for the matrix B in question 3 and 
hence reduce B to diagonal form. 
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Rectangular Matrices 


The treatment of matrices so far has been confined to (2x2) and 
(3 x 3) matrices and the use of the column matrix to represent a 
mathematical quantity, i.e. a vector. When the number of columns 
and rows are equal, the matrix is said to be square and the 
treatment of square matrices forms the greater proportion of the 
subject at this level, but other matrices, with the number of columns 
not equal to the number of rows, called rectangular matrices, 
must receive some consideration. They often arise in the treatment 
of linear equations. 


y = a l x l + b l x 2 + c l x 3 
/ = a 2 x l + b 2 X 2 +c 2 x 3 

In matrix form this becomes 


y 


b , c, 



y. 


-«2 b 2 c 2 . 


*2 




-*3- 


The matrix 


b i 


has three columns but only two 


l a 2 b 2 c 2 j 

rows, it is called a (2x3) rectangular matrix. (Note: Rows are 
always given first, columns second.) 

The product of the (2 x 3) matrix with the (3x1) column matrix 
yields a (2 x 1) column matrix. Inspection shows that two matrices 
can be added only if they have the same number of rows and 
columns, i.e. they are said to be conformable for addition. 

Two matrices can be multiplied together only if the number of 
columns of the first is the same as the number of rows of the 
second, and again they are said to be conformable for multi¬ 


plication. 

For matrices which are not square (i.e. rectangular) it is im 
possible to define a determinant or an inverse. 
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Introduction to Matrices 
Example ( i ) 



1 

3 

-1 

2' 


i 

1 

0 

2 

Add A = 

3 

2 

2 

0 

to B = 

2 

1 

0 

3 


.1 

3 

5 

-2. 


.0 

4 

-1 

3. 


Example (ii) 
Multiply C = 


A + B = 


1 4 2 
-13 1 


4-1 4 

3 2 3 

7 4 1. 


by D = 


2 -I 
4 3 
2 5j 


CD = 




'22 

21' 


.12 

15. 


.(-2+12 + 2) (1+9 + 5) 


C is a (2 x 3) matrix, D is (3 x 2) matrix, and the product is a 
(2 x 2) matrix. 


Examples XX111 

\2 3-1 2 
.3 4 6 1 


A = 


B = 


6 3 
2 1 


C = 


4 0 
-I 6 


2 1 
3 8 
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'2 

3' 


‘ 2 

1 

0 

3" 


"2 

1 

2" 


4 

1 


3 

2 

1 

6 


3 

0 

6 

D = 

6 

8 

£ = 

-1 

2 

0 

0 

F = 

2 

1 

8 


.9 

2. 


4 

3 

1 

2. 


.4 

1 

6. 


G = 


1 -1 


1 Describe A, B. C, D. £, F and G as (nxm) matrices, (//infr 
/4 is a 2 x 4 matrix.) 

2 State which matrices are conformable for addition and give 
the sums. 


3 State which matrices are conformable for multiplication and 
give the products and class of matrix (n x m) involved. 

4 Expand 


x' 


-1 

o 

UJ 

_ 


"x 



.0 1 4. 






.1. 


and state the kind of geometrical transformation it represents. 
5 What kind of transformations are represented by 


V 


1 0 3' 

X 



o 

1 

1 

4^ 

y 

_i. 


(2 x 2) Matrices and Quadratic Forms 

The transpose of a column matrix is a row matrix. 


If A = 


then A' = [x y] 


y J 


A 1 .A 


[X 



= x 2 +y 2 
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Introduction to Matrices 


Next we examine the matrix product which uses this property: 


[X y] 


5 

2 


yl 


= [x >•] 


(5x+2y) 
,(2x+3y)j 
= 5x 2 + (2xy+2xy) + 3y 2 
= 5x 2 + 4xy+3y 2 


The equation 


5x 2 +4xy + 3y 2 = 5 
can be written in matrix form 

[5 21 


[x >-] 


2 3 


y 


= 5 


The general homogeneous quadratic form ax 2 + 2hxy + by 2 = 1 
can therefore be expressed in matrix form 


[x y] 


a W 

X 

.b b_ 

.y. 


= 1 


The equation ax 2 + 2hxy + by 2 = 1 represents a conic section and 
by choosing suitable values for the elements a, b, h of the matrix, 
the matrix form can be used to represent the circle, ellipse, hyper¬ 
bola and so on. By a suitable transformation the axes of the conic 
can be rotated if this is necessary. 

The general quadratic form 


lx y] 


can be written more briefly X'AX since the row matrix is the 


a h 

X 

h h . 

. y. 


transpose of the column matrix and A = 


If A = 

LO 

the equation 


0 

b 2 


[x y] 


Ci ' 

N> 

o 

X 

.0 b 2 . 

. y . 


= a 2 b 2 
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becomes 


or 


a 2 x 2 + b 2 y 2 = a 2 b 2 
x 2 y 2 

—x -}- = 1 i.e. the ellipse. 

b or 


[x y] 

a h 

> 


h h. 



the equation 
= 1 represents the equation of the tangent to 


If the column matrix 

x‘ 

is replaced by 

x 1 


.y. 




the conic on page 110 at the point (x\ /) or 


X'AX' = 1 where X' = 


(3 x 3) Matrices and Quadratic Forms 


Another matrix form using a column matrix 


x 


and its transpose is 


y 

1 



'a h g~ 

X 

[x y 1] 

h b f 

y 


bo 

.1. 


The (3x3) matrix is observed to be symmetric since the 
elements in corresponding places on either side of the leading 
diagonal are equal, i.e. its transpose is equal to itself ( A 1 = A). 
Multiplication of the three matrices gives 


[x y U 


ax by g 
bx by f 
gx fy c 


= 0 


-* ax 2 + by 2 + 2bxy + 2gx + 2fy+c = 0 

which is the general equation of the second degree and represents 
a conic. 


Ill 



















































Introduction to Matrices 
The special case 


gives 



'1 

0 

z 

x" 

[x y 1] 

0 

1 

f 

y 


-g 

/ 

c_ 

.1. 


' x+Z 

[x y 1] y+f =0 
_gx +fy + c. 

-* x 2 +gx + y 2 +fy+gx+fy + c = 0 


i.e. 


x 2 +y 2 + 2gx + 2fy+c = 0 

familiar to most as the general equation of the circle. The tangent 
to the circle at the point (x'y') is 



'1 

0 

8 

V 

[x y 1] 

0 

1 

f 

y' 


.g 

/ 

c_ 

.1 _ 


xx' + yy'+g(x+x')+/(y+y')+c = 0 


If X = 

X 

, then X x = [x y z] and A = 

a h g~ 


y 


h b f 


_ z_ 


..g f c. 


then 


X'AX = 1 gives an homogeneous equation of the second degree 
in x, y and z representing a quadratic surface. 

X'AX = 1 on multiplication gives 

ax 2 + by 2 + cz 2 + 2fyz + 2gzx + 2hxy = 1 


Examples XXIV 

1 Find the quadratic form equivalent to the matrix product 


[x 


y] 


4 3 
3 2 


x 

.y. 
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2 Find the quadratic form equivalent to the matrix product 


3 

4 

5 

6 
7 


[X y] 

1 2" 

x' 


.2 1. 

.y. 


Express the quadratic form x 2 + 2xy+y 2 in matrix notation. 
Express 2x 2 + 6xy + 6y 2 in matrix notation. 

Express 2x 2 + 3xy + 6y 2 in matrix notation. 

Express x 2 + y 2 in matrix notation. 

Express x 2 + y 2 = 25 in matrix notation. 


8 Find the quadratic form equivalent to the matrix product 


[X y] 

2 O' 

x' 


.0 5. 

.y. 


9 Solve the matrix equation [x 2][ 2 0 

■ • U J 

10 Find the quadratic form represented by 


= 0 . 



'1 

0 

3‘ 

x" 

[x y 1] 

0 

1 

2 

y 


.3 

2 

-3. 

.1. 


= 0 . 


11 Give the circle x 2 +y 2 + 6x-8y + 9 = 0 in matrix form. 

12 Give the equation of the tangent to the circle 

x 2 +y 2 + 6x — 8y+9 = 0 
at the point (x\ /) in matrix form. 

give the quadratic form of the 


13 If A = 

'1 0 3' 

, x = 

X 


0 1 2 


y 


.3 2 -3. 


_ 2 _ 


matrix equation X'AX = 0. 
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Introduction to Matrices 

General Definition of a Matrix 

Throughout all the work done on (2 x 2), (3 x 3) and rectangular 
matrices, the position of the element of a matrix is of the greatest 
importance. The position of the element is completely defined by 
its row and column position. 

Definition 

A rectangular array of numbers (called elements) containing in 
rows and n columns is called a matrix. 

The elements are enclosed in square brackets (used throughout 
this book) or in curved brackets (widely used in many modern 
books). Each element is given two suffixes, the first to denote the row 
and the second the column. This convention is always followed. 

a i 1 a l2 (, I3 ■ a ln 
a 2i a 22 a 23 ■ ■ ■ a 2n 


_ l tl m 2 O m3 * * • 

In the above arrays a 23 means the element in the second row 
and in the third column; a mn means the element in the /nth row 
and the nth column. The general element is usually denoted by 
a„ or a u . 

Abbreviated Notation 

The notation for a matrix can be abbreviated still further and for 
theoretical work the contraction makes for simplicity of expression. 
The matrix shown in the previous paragraph can be denoted by 
K'll or [ay] or ( a ij) or A. For some purposes it is useful to speak 
of one of the terms of the matrix and here another notation is 
used. If [A] or (A) is the matrix A; then [A] 0 or (A) (j means the 
element in the ith row and j'th column of the matrix A. 

The sum of two (m x n) matrices A and B is obtained by adding 
corresponding elements as 

A = IKH B = || b,j\\ 

a +b = KII+IIM = K-+M 

114 


or 


An 

°12 

<* 13 • 

• * ^1 n 


°22 

<*23 * 

■■ «2» 

a 3i 





\ ^m2 ^m3 • • • tl mn j 
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But since the sum (a,j+b,j) is a sum of numbers these numbers obey 
the usual commutative and associative laws for numbers , we can 
write 

A + B = || fly- +b,j || = ll&y+flyl = B + A 
{A + B) + C = IK-+M + IM = ||K+*>y)+ c ull 

= |fly + (6y + Cy)|| 

= ll a ul + ll b u +c ul 

= A+(B + C) 

If k is a scalar then 

kA+kB= \\kaijW + ||/cbj = \ka u +kb u \ 

= |*K+MI 

= k(A + B ) 

Multiplication of Matrices 

c = A.B 

First we remember that matrices A and B must be conformable, 
i.e. the number of columns of A must be equal to the number of 
rows of B. 

With A and B both (2x2) matrices each element of C, their 
product, was composed of the addition of two products. 

With A and B both (3x3) matrices the elements of C were the 
result of addition of three products. 

If A is an (m x n) matrix and B is an (n x I) matrix they are con¬ 
formable and C will be a (m x /) matrix and each element of C 
will be the sum of n products. 

k -n 

[AB] (j = Y. Cl ik-bkj 

k ts l 


or 

k = n 

Cij = X a >i‘ bkj = ^ n bij+a n b 2 j+a i3 b 3 j . a in b nj 

k= 1 
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Summary of the Laws of Matrix Algebra 

The laws of addition and multiplication given on page 56 for 
(2 x 2) matrices can now be extended to all matrices which are 
conformable and remembering that, in general, multiplication of 
matrices is nor commutative. Some additional relationships 
covering inverses and transposes must be given. 

The following laws are restated for all conformable matrices. 

(i) A + B = B + A (Commutative Law for addition.) 

(ii) (A + B) + C = A+(B + C) (Associative Law for addition.) 

(iii) A + 0 = O + A = A (Null or Zeromatrix is the Identity 

element for addition.) 

(Associative Law for multiplication.) 
(Distributive Law for multiplica¬ 
tion.) 


(iv) (AB)C = A(BC) 

(v) (A + B)C = AC+BC 
and 

C(A + B ) = CA + CB 


(vi) A. I = IA = A 


(Note: the order of multiplication 
must be adhered to since multiplica¬ 
tion is not commutative.) 

(Unit matrix is the Identity element 
for multiplication.) 

The inverse of matrix A is denoted by A -1 but the Transpose of 
A can be denoted A' or by A' either notation being acceptable. 
Since A' or A' results from the interchange of rows and columns 

(A 1 )' = A 

and if the product AB exists then it was shown earlier that 
(AB)‘ = B'A' 
hence if ABC exists then 


(abcy = (ab. cy 

= C.(AB)' 

= C.B'.A' 

It can be shown also that since A.A~ 1 = I 
(A-'y = (A ')~ 1 

and if the product AB exists and A and B are non-singular then 
(ABC 1 = B-'A- 1 
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Examples XXV 


A = 

'2 3 

B = 

-2 2 

c = 

2 -f 


.4 7. 


-3 2 . 


.3 4. 


1 Find A\ B‘, (AB) 1 , B 1 .A 1 . 

2 Find C', (ABC)', C.(ABy, C.B'.A'. 

3 Find AB, (ABV 1 , B~ 1 A~ 1 . 

4 Find (C~ *)', (C') -1 and state your conclusion. 


5 P = 

'1 

0 

3' 

Q = 

'3 

2 

r 


2 

1 

4 


4 

2 

2 


.0 

1 

3. 


.1 

3 

0 . 


(a) Find P‘,Q‘,(PQ)', Q'.P'. 

(b) Find P -1 , Q- 1 , (PQ)' 1 , 

6 Since A.A -1 = / and /' = I, prove that (A -1 )' = (A') -1 . 

7 Why is (A — B)(A + B), in general, not the same as A 2 — B 2 ? 

8 Expand (I — A)(I + A). 

Linear Dependence and Rank 

The concepts of Rank and Linear Dependence are of the utmost 
importance in the further applications of matrix algebra. 

Linearly Dependent Quantities 

A set of “quantities” x„ x 2 , x 3 .. .x„ is said to be linearly dependent 
if there is a linear combination of them equal to zero, 

a,.*, +02X2 + 03*3+ • • a « x n = 0 

with coefficients a„ a 2 , a 3 ...a„ in the field of numbers in which 
one is working and not all these coefficients being zero. 

If the only relation possible is given when a t = a 2 = a 3 ... a n = 0 
then the quantities are said to be linearly independent. 

Example (i) 

The quantities x + 2y and 3x + 6y are linearly dependent because 
-3(x + 2y) + (3x + 6y) = 0 
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Example (ii) 

The quantities ( —x4-2y4-3z), (x-4y-13z) and (-3x + 5y+4z) 
are linearly dependent since 

7( —x4-2y4-3z)4-(x —4y—13z)—2( —3x + 5y4-4z) = 0 


Example (in) 

The quantities (-x + 2y-3z), (2x-y + 4z) and (3x + 4y + z) are 
linearly independent because if 

a( — x4-2y — 3z)4-6(2x — y + 4z) + c(3x + 4y + z) = 0 

this leads to the three relations 

x( — a + 2b + 3c) = 0 => — a + 2b 4- 3c = 0 

y(2a — h + 4c) = 0 =♦ 2a — b + 4c = 0 

z( —3n + 4h + c) = 0 => — 3a4-464-c = 0 

Eliminating a gives 

26 + 8c = 0 


and 


36 4- 10c* = 0 

which is only possible if b = 0 and c = 0 and therefore a = 0. 

Hence by definition (-x + 2y-3z), (2x-y + 4z), (3x4-4y4-z) 
are linearly independent. 


Example ( iv ) 


The three vectors 

'3 

» 

2' 

> 

4 


2 


0 


2 


4 


2 


.3 


are linearly independent because 


the relation a 

'3' 

+ 6 

'2 

4-c 

4" 

= 

O' 


2 


0 


2 


0 


.4. 


2_ 


.3. 


.0. 


"3 a~ 


'2b' 


’4c" 


'O' 

2 a 

4- 

0 

4- 

2c 

= 

0 

4 a 


2b 


3c. 


0 


leads to 
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'3a+ 26+4c' 


'O' 

2a +2c 

= 

0 

_4a + 2b +3c. 


.0. 


for these two column matrices to be equal their corresponding 
elements must be equal 

3a + 26+4c* = 0 
2 a 4-2c = 0 
4a4-264-3c = 0 

and the only possible solution is 

a = b = c = 0 


Examples XXVI 

1 Show that (x 4- 2y) and (3x-4y) are linearly independent. 

2 Show that the quantities (x4-3y) and (2x-6y) are linearly 
independent. 

linearly dependent? 


3 Are the vectors 

2 

and 

T 


4 


.3. 




4 Are the quantities (2x — 3y + 4z), (3x + 4 y—z) and (— 6x — 5.y — 2z) 
linearly dependent? 

5 What is the condition that (ax + by) and ( cx+dy ) are linearly 
dependent quantities? 

6 If p = (2x+y + zX q = (3x-4y + 6z) and r = (4x-9y+ll, 

find the relation between p, q, r for linear dependence. Calculate 
2 1 1 and state your conclusion. 

3-4 6 

4 -9 11 

7 If n = (2x+y + z), in = (x + 2>- + 2z) and p = (x + 3y+2z) are 
the quantities n, m and p linearly dependent? 

linearly dependent? 


8 Are the vectors 

"2' 

> 

'-1' 

* 

"3' 


1 


2 


4 


.4. 


.-7 


. 1 . 
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Introduction to Matrices 
9 Calculate the value of 


2 — 1 3 
1 2 4 


and compare with the result 


4 -1 1 


of question 8. State your conclusion. 


£ 


Rank of a Matrix 

Definition 

The hank of a matrix A is defined as the maximum number of 
linearly independent rows of A. 


It can be shown that rank defined in this way is also equal to 
the order of the largest non-zero (or non-vanishing) minor of the 
determinant of that matrix. 


Matrix A = 
The rows (3 


3 2 . 

1 4, 

2) and (1 4) are linearly dependent if 

a(3 2)+b(l 4) = 0 


*■* 3 a+b = 0 a = 0 
-* 2a+4b = 0 b = 0 


The rows are linearly independent, hence the rank is 2. 


Matrix B = 


a - 2, b = 1. 


3 1 
6 2 


and a(3 1) + b{6 2) = 0 is satisfied by 


The rows are linearly dependent and the rank is less than 2 
but equal to 1 because there are non-zero elements in the minors. 
Only the null matrix has rank of zero. 

Notice that det|/l| = 10 hence rank is 2 but det|B| = 0 hence 
rank is less than 2, i.e. 1. 

Consider the rows of the matrix P = T3 4 1 " 

4 3 2 

2 1 4. 


) 
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a(2 4 l)+b(4 3 2) + c(2 1 4) = 0 


3a + 4b + 2c = 0 
4a + 3b+c = 0 
a + 2b + 4c = 0 


Section IV 


From the first two equations —5a —2b = 0 

From the last two equations — 15a — 10b = 0 

and therefore the only possible solution is a = b = c = 0. The 
three rows are linearly independent, the rank is 3. 

Had we used columns instead of rows for testing we would have 
found 


resulted in 


'3" 


4 


T 


0" 

4 

+ b 

3 

+ C 

2 

= 

0 

2. 


_ 1 _ 


.4. 


.0. 


3a+4b+c = 0 
4a + 3b+2c = 0 
2a + b+4c = 0 


and again the only solution is a = b = c = 0. 


The determinant 


3 

4 
2 


4 

3 

1 


1 

2 

4 


— 20 and it is from a (3 x 3) matrix. 


Therefore it is of order 3, hence the rank is 3. 


The matrix Q = 


3 5 -15 

4 3 2 


has its three rows linearly depen- 


J 1 4. 

dent because they can be expressed in the linear relationship 


a(3 5 —15)+/>(4 3 2)+c(2 1 4) = 0 

1 . v • 

where a = 1, b = — 3^ and c = 5\. The rank is less than 3 but 
since any two of the rows are linearly independent then the rank 
is 2. 

Notice that the determinant of Q is 0 but any (2x2) sub-matrix 
has a non-zero determinant, the rank is 2. 


E 
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Introduction to Matrices 


The 3x4 rectangular matrix R 


3 

4 
1 


2 1 
3 1 
2 1 


has only 3 rows 


but 4 columns. 

The four column vectors are linearly dependent since if 


'3' 


"2' 


T 


'8' 

4 

+ b 

3 

+ c 

l 

+ d 

2 

.1. 


2 


l 


6 


then a = 1, b = — 7, c = 19 and d = —1 is a possible solution 
and therefore by definition the four columns are linearly dependent, 
however any three columns are linearly independent. The three 
row vectors [3 2 1 8], [4 3 1 2], [1 2 1 6] are, how¬ 

ever, linearly independent. Hence the rank calculated by columns 
is 3 and the rank calculated from rows is 3. This illustrates an 
important property of any matrix that 


Row Rank = Column Rank 

The matrix 


S = 



2 1 
0 4 
6 8 
-2 4 


3" 

1 

2 

1 


is a square (4 x 4) matrix of order 4. 

Hence its rank could be 4, 3, 2, or 1. 

Calculation shows that the four row vectors are linearly inde¬ 
pendent, the rank is 4. 

The diagonal matrix 


D = 


0 0 0 
b 0 0 
0 c 0 
0 0 d 
0 0 0 


0 

0 

0 

0 

e 


122 


Section IV 


has one non-zero element in each row on the leading diagonal 
and inspection shows that the rank must be equal to the number 
of non-zero elements along the diagonal. Rank of any diagonal 
matrix can be seen by inspection. 

It is a property of a matrix that the Rank remains invariant 
during any of the elementary transformations previously 
defined. By reducing a large matrix using elementary transforma¬ 
tions to an equivalent diagonal form, its rank can be seen by 
inspection and is equal to the number of non-zero elements on 
the leading diagonal. 


Example 

The method will be illustrated by using the (3 x 3) matrix 
M=T3 2 2"| and the transformations will be listed. The 

1 1 4 

2 3 2. 


sign ~ means equivalent to. 



'3 

2 

2 " 


1 

1 

4 


1 

3 

2 . 


'1 

§ 

r 


1 

1 

4 


.1 

1 

1. 


"1 

f 

2 

7 


0 

i 

H 


.0 

s 

7 

i 

7 


'1 

0 

— I 

•"'w 

0 

-8 

\ 


.0 

s 

5 



"1 

0 

— | 


0 

-8 

\ 


.0 

0 



Divide row 1 by 3 and row 3 by 2. 


Subtract row 1 from row 2 and row 3. 


Subtract 2 x row 2 from row 1 and 
10x2 row 3 from row 2. 


Add jg x row 2 to row 3. 


Add 18 x row 3 to row 1. 
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Introduction to Matrices 

Ti o o 
0-8 0 
.0 0 f 


This matrix is now in diagonal form, its 
rank is 3 and since the"rank has remained 
invariant throughout, the rank of matrix 
M is 3. 


Examples XXVII 

\ What is the rank of A = 

2 What is the rank of B = 


2 -1 
3 4 


3 What is the rank of C = 


4 What is the rank of D = 


5 What is the rank of £ = 


6 Reduce S = 


1 

2 

1 


2 

0 

6 


3 2 4 
2 0 2 
L4 2 3 

2 3 

3 5 
2 4 
2 4 

-12 3' 

1 -4 -13 
-3 5 5J 

1 3 by elementary transformations 

4 1 

8 2 


and find |D|? 

and find |£|. 


-1 -2 4 1. 
to an equivalent diagonal form and deduce the rank of S. 


Rank and Consistency of Linear Equations 

The concept of rank can be used to state the condition for con¬ 
sistency of a group of simultaneous linear equations, m linear 
equations in n unknowns are consistent when, and only when, the 
rank of the matrix of the coefficients is equal to the rank of the 
augmented matrix. The Augmented Matrix is the matrix of co¬ 
efficients to which the column of values (18 10 29 in the 

example below) has been inserted as a final column. 
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Example 


4x+3 y+z = 18 
2x+y + 3z = 10 
5x + 7y-2z = 29 


The matrix of the coefficients is 


4 

2 

L5 


The rank of this matrix is 3. 


The augmented matrix is 


3 1 
1 3 
7 -2 


7 -2. 

18 ‘ 

10 

29 


and the rank of this 


matrix is 3. 

Since the rank of the augmented matrix is equal to the rank of 
the coefficient matrix the set of equations is consistent, i.e. these 
equations are satisfied by some set of values of x, y and z. 

But if 

4x + 3y + z = 0 
2x+y+3z = 0 
5x + 7y —2z = 0 


the matrix of coefficients is again 


'4 3 
2 1 


.5 7 -2. 

The above set of homogeneous equations can be expressed as 
AX = 0 where A is the matrix of coefficients and X is the column 


matrix 


z_ 

Since \A 0 then X = 0, i.e. there is only the Trivial solution 
x = 0, y = 0, z = 0. 

The application of matrix algebra to the study of Linear Equa¬ 
tions is an important branch of Modern Algebra and is dealt with 
more fully in a number of the books in the bibliography (6, 7, 9, 
10, 15). 
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Rank and Image Space 

A line is one-dimensional, an area is two-dimensional, a volume 
is three-dimensional, a space of n dimensions is spanned by hyper- 
volume ! 

Three linearly independent vectors are needed to span three- 
dimensional space, and two linearly independent vectors to span 
two-dimensional space. 

From this we can give some geometrical illustration of rank. 

The (2 x 2) matrix represents a transformation and if the rank 
is 2 the image produced will be in two dimensions. If the rank is 1, 
then the image will be one-dimensional, i.e. a line not an area. 
This was shown on page 36 in diagram 16. 

The (3 x 3) matrix represents a transformation in three dimen¬ 
sions—generally the unit cube transforms to a parallelepiped, if 
the rank of the matrix is 3. If the rank is 2, the image of the unit 
cube is in two dimensions only. (This follows because if the rank 
is only 2 the determinant of the matrix is zero hence the volume 
of the parallelepiped is zero also, i.e. it has become an area.) 

Generally then the rank of a matrix is the dimension of the 
image space. 


History of Matrices and Determinants 

The form in which we use matrices is due largely to the work of 
Cayley, 1841, and Sylvester, 1839, in the early part of last century, 
but the use of determinants with a poorer form of symbolism had 
a much longer history mainly in connection with the solution of 
simultaneous equations. 

Eight hundred years ago the Chinese were solving linear equa¬ 
tions by using rods on a calculating board and developed the idea 
of elimination which was virtually a method of subtracting rows 
and columns of a determinant. 

A Japanese mathematician, Seki Kowa, in a book written 
about 1683, showed some understanding of determinants which 
he used for elimination of unknowns from sets of linear simultane¬ 
ous equations. 

The theory of determinants was first developed in Europe by 
Leibnitz again in relation to the solution of simultaneous equa¬ 
tions but the symbolism was unlike our own. Vandermonde 
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in 1771 laid the foundations of modern theories of determinants 
as functions independent of any applications to linear equa¬ 
tions. 

Laplace (1772) showed how to expand a determinant in terms 
of its minors and Lagrange in 1773 treated determinants of the 
second and third orders. Gauss (1801) advanced the theory of the 
subject and in 1812 Binet showed how to multiply two matrices. 
About the same period Cauchy finally established the use of the 
word ‘determinant’ in its present sense and made great contribu¬ 
tions to the advancement of the theoretical treatment of deter¬ 
minants. The work and contributions of the great mathematician, 
Jacobi (1839), in this field are remembered by naming certain 
forms of matrices after him. From this time (1839) forward the 
symbolism and theory of matrices received its greatest contribu¬ 
tions from the work of Cayley (1821-1895) and Sylvester (1814— 
1897). 

In 1850, Sylvester used the word matrix for arrays of numbers 
from which determinants could be formed and the theory of 
determinants is now regarded as a branch of matrix theory. 
During the early part of the nineteenth century determinant 
theory was applied intensively to the study of systems of linear 
equations, but the concepts of rank and linear dependence of a 
matrix which Sylvester developed, have produced a new outlook 
in the treatment of linear equations; the rapid development of the 
high speed computer from 1946 onwards and the wide applica¬ 
tions of linear programming (1947) have hastened the change in 
approach. 

In 1841 Cayley introduced the determinant notation which is 
now in general use and has been since he first announced it, then 
in 1858 he defined the square matrix in the form which has become 
standard and he continued to contribute important theory of the 
subject which the student will meet later. 

Sir W. R. Hamilton discovered the principle of Quaternions in 
1843 and published the full treatment in 1853 and it is in this 
connection that he is usually remembered but he made extensive 
and valuable contributions to matrix theory which he used in his 
own work. 

Calculations with matrices were much simplified by the use of 
elementary transformations due to Grassman, 1862, and Kronecker, 
1866, while Frobenius (1870-1911) developed ideas in equivalence 
in the theory of matrices. 
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Introduction to Matrices 

So far in this Introduction to Matrices , the elements have been 
in the field of real numbers but elements in the complex number 
field can be used and the subject pushed deeper: complex matrices 
can be used in the representation of quaternions. 
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Answers 

Examples I 


1 (i) 

y 1 ' 


'2 

3' 


(H) 

>i‘ 


10 12 ' 

r*r 


. y 2 . 


.4 

2. 

.*2. 


.yi. 


.25 16. 

U 2 . 

(iii) 

x' 


'3 

2' 

’ X 

(iv) 

x' 


4 “ 3 1 

X 


y 1 


2 3 


y\ 


y 


2 -1 


(i) x' = 2 x+y 
y' = 4x + 6y 


(ii) x' = 3 x+y 
y' = 2 y 


3 C = 

6 

9 

D = 

'2 3' 


.6 

11. 


.4 3. 


4 B + A = 


5 P + Q = 


(P+Q)+R = 


6 9 

6 11 . 

3 5' 

4 10. 
"3 

4 


-4 + B = B+A 


Q + R = 


P+(Q + P) = 


5 

10 

3 


4 7 

7 11J 
4 7 

7 11 


6 q+p = 

p+q = 

7 (p -0 = 


3 5 

4 10 
2 4 

6 5 
1 1 
-2 2 


(P-Q)-P = 



Q + P = 

P+Q 


P + Q = 

Q + P 


(Q-P) = 

0 O' 



.0 3. 

-1 


\2 


-5 1 


P-(Q-R) = 


3 

1 3 


y 
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(P-Q)+R = 


2 3 
1 3 


Examples II 
1 i4.B = Til 16 
.19 26 J 
3 -6 
-1 2 
11 16 
.19 26 
’ 7 4 

.15 10J 
(AB)P = A(BP) 

4 B./4 = [ 4 15 

.10 33 
'll 7 
.19 8 
(B/4)P= B(/4P) 


2 P.Q = 


3 A.B = 


B.P = 


A.P = 


5 /4 + B = 


P4 = 


3 5 

4 10 

4 15 

7 15 J 


PA + PB = 


B + P = 


11 25 
13 25 
2 4 


4B = 


6 5 
11 16 
19 26 


B.A = [ 4 

.10 

fi.P- 

(/l.B)P = 

-4(B.P) = 


(BA)P = 

B.(/4P) = 


P(A + B) 


PB = 


1 
.6 

P(A + B) ■■ 


A(B + P) = 

r 11 


AP = 


19 


15' 

33. 

5 ° 1 . 

10 0 . 

'59 38' 

.97 64 

'59 38' 

97 64. 

49 23' 
.109 53. 
49 23' 
*[l 09 53. 

11 25 
IL13 25. 
10 ' 

10 . 

= PA + PB 

'22 23' 
’ .38 34 
7] 

8 . 
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AB+AP= 


Introduction to Matrices 
22 23' 

.38 34. 

A(B + P) = AB + AP 

6 Matrices do not obey commutative laws of multiplication. 
AB # BA ( AB)P # (BA)P 

Examples III 


1 2/4 = 

4 6 

5/4 = 

10 15 

-34 = 

‘ -6 -9' 


.2 12. 


. 5 30. 


OO 

7 

r*"> 

1 


kA = 
4 B = 


2k 3 k 
k 6k 
4 8 

12 16 


kB = 


k 2k 

L3fc 4* 


-2B = 


-2 -4 
-6 -8 


3 2P = 

'4 —2 

nP = 

2/j -n 

— 3P = 

-6 3 


.6 -6. 


,3n —3n. 


. -9 9. 


4 A + A = 


4 6 

2 12 


= 2/4 


Examples IV 


2 DK = 


3 MB = 


IM = 


6 8 
2 12 

2 0 

-1 0 

2 0 

-1 0 


= /CD 


KE = 


BM = 


IB = 


2a 2/> 
2c 2d. 
2 0 
-1 0 
1 0 
1 3 


= EK 


4 4-/ = 

1 3' 

A —31 = 

-1 3' 


.1 5. 


1 3 


2 B-I = 


1 4 
6 7 


B-/c/ = 


(I-*) 2 

3 (4-k) 
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/4/j.VM-cr.v 



15 

10 


Examples V 



X 

= 

' 2 

3' 

V 


.y . 


3 

5. 

- V- 


2 

3 



3 0 
0 1 
-1 
0 


x 


This represents a shear. 


(2.3) becomes (-2,3) 
(3. 5) becomes (- 3. 5) 


5 The axes are rotated clockwise through —60°. 

6 

7 All lengths are doubled. 

8 Both are shears, but the second is 3 times larger in the x 
axis than the first. 


9 


Examples VI 

1 -2 

2 32 

3 0 singular 


4 A"' = 


5 -7 
-2 3 


-2 0 
0 -2 


20 4 

10 22 


AA ~ 1 = A~'A = I 


2 3 
4 5 


" = -i 


5 -3 
-4 2 


6 

'll 4 


67 

41 

P~' = 

1 

n 


1 8 


.29 

19. 


.-3 5. 


TERM 
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Introduction to Matrices 


8 


13 


3 -4 

-2 7 


1 


(ps-qr) 


10 -1 


8 


s -q 
-r pi 
5 


(ps-qr) 0 


11 


.-3 -2. 

P'x' = 13, y' = 9 
Q'x' = 33, / = 23 


P"x" = 2, / = 1 

Q"x" = 4, y" = 3 

Conclusion: The images of F and Q' due to the inverse matrix 
are the original points P and Q. 

12 A non-singular matrix represents a reversible mapping 
and the Inverse of the matrix reverses the mapping. 

Examples VII 


1 

x' 

= 

4 5 

X 




.2 3. 

.y. 


0' is (0, 0) P' is (4, 2) Q' is (9, 5) R' is (5, 3) 
The new area is 2 sq units. 

4 5 

A = =2 New area = & x Old area 

2 3 

Area A'B'C'D' = 5. Area ABCD 


3 2 
2 3 


= 5 


3 ABC' 
ABC 

2 3 

3 2 


3 2 
2 3 
3 2 
2 3 


Examples VIII 


1 

'2 

3 

x' 

= 

18' 


.3 

5 . 



.29. 
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Answers 


2 2x-3 y = 4 
2x + 3y = 6 


3 P.X = 

2x + 3y 

P~ l (PX) = 

X 


.3x + 5y. 


.y. 


4 

5 

6 
7 


x = 3 


x = ^ = 67 


C = 


Zl = 


11 1 

.15 2. 
2y, + 3y 2 


z 2 = 3y t + 4y 2 


y = 4 
y - i 

y = = 74 

y, = x, + 2x 2 z 2 = llx,+x 2 

y 2 = 3x 1 — x 2 z 2 = 15x, + 2x 2 


Examples IX 


1 

2 


P' = 

P' = 


2 —2 
5 4. 

2 31 


3 5 


3 

e _1 

= 

‘ cos a - 

sin a 

Q‘ = 

' cos a - 

sin a ’ 




sin a 

cos a. 


sin a 

cos a. 

4 

R -1 

= 

cos X 

sin x n 

R' = 

cosx 

sin x' 




— sinx 

cos x^ 


— sin x 

cosx. 


5 Q and R are orthogonal. 


6 A' = 

3 2 

B' = 

2 4' 

(AB)' = 

'22 24 1 


A 5. 


.3 7. 


.37 41. 


A'B' = 

12 

26' 

B'A' = 

'22 

24' 


.23 

51. 


.37 

41. 


conclusion (AB)' = B'A' 


7 A = 7 
B = 2 


A' = 1 
B' = 2 
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Introduction to Matrices 


8 A' — B' = 

'5 6 

1 

s 

II 

'5 6 


.7 12. 


.7 12. 


Examples X 


-5 5 


= 4 + / 


4 1 
-1 4 

3 (2—3/)(4 — 5i) = 


\_ 

13 




6 r 


-5 12 

-12 -5 
0 13 
-13 0 

cos 0 sin 0 
— sin 0 cosO 


-7 -22 

22 -7 

= -5+12/ 


= -7-22/ 


1 2 + 3/ = 

' 2 3" 

3 + 2/ = 

-3 2' 


.-3 2. 


1 

1 

u> 

(2+ 3/)+ (3 + 2/) = 

r 5 5] = 5 + 5/ 



7 2 , = r, 

cos 0, 

sin 0, 

; 2 2 = r 2 


. - sin 0, 

cos 0, 



cos 0 2 sin 0 2 
— sin0 2 cos 0 2 


z,~ 2 = r,r 


1'2 


cos (0,+0 2 ) sin (0,+0 2 ) 
-sin (0,+0 2 ) cos(0,+0 2 ). 


Examples XI 


3 OC 


4 

3 


4 OPRQ is a Parallelogram 

5 ORAQ is a Parallelogram 
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Answers 


Examples XII 


1 AH = 

2 KB = 

3 LC = 

4 DN = 



HA = 


BK = 


CL = 


N£> = 


-2 1 
3 4 
4 15 

2 4 

9 -2 

-5 3 

10 2 
9 -18 


5 

6 


T 1 = 

0 

0 

T A = 

0 

9 

AT = 

0 

9 



.0 

0 . 


.0 

-6. 


.0 

-6 

I 

F 2 = 

0 

4 ' 

FB = 

4 

-4 ‘ 

BF = 

0 

14' 



.0 

4 . 


.4 

-4 . 


.0 

0 . 



II 

' —2a 

-2b' 

PA = 

-2a 

-2b 

AB = 

O 

1 

ON 


. —2c 

-2d. 


. -2c 

-2d. 


. -6 0. 


BA = 

0 -6' 

BP = 

3c 3d' 

PB = 

3b 3a 


. -6 0. 


.3 a 30. 


i- 

sx 


Examples XIII 


>/' 


'2 

3 

4' 

"*r 

)’2 

= 

1 

-2 

3 

^2 

3’3- 


_3 

2 

-2. 



2 = 3x, + 2x 2 + x 3 

y 2 = — 2x, +3 x 2 + 6x 3 
y 3 = x, — 2x 2 + x 3 
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Introduction to Matrices 


3 P+Q = 

-1 

-J 

OJ K) 

_1 

e+p = 

‘4 4 2' 

7 4 3 


.3 6 3. 


.3 6 3. 


Q-P = 


2 0 0 
1 0 1 
-1 0 -lj 


4 3 B = 

"3 

6 

3' 

kB = 

' k 

2k 

k' 


9 

6 

9 


3k 

2k 

3k 


.6 

3 

6. 


.2k 

k 

2k. 


5 P+Q = 

"2 3 2' 

(P + 0 + P = 

"3 3 3' 


3 1 3 


5 2 6 


.6 3 5. 


00 

r- 

Q + R = 

'2 12' 

P+(Q + R) = 

'3 3 3‘ 


4 2 5 


5 2 6 


.4 4 6. 


1 

—1 

LA 

OO 

L_ 


i.e. (P + Q) + R = P+(() + /?), this verifies the Associative Law 
of Matrix Addition. 


Examples XIV 

1 P.Q = 


2 B.l = 


"12 

9 

6' 

Q-P = 

"11 

13 

7* 

18 

13 

8 


14 

18 

10 

.20 

16 

10. 


.12 

11 

6. 


B + kl = 


"1 

1 2" 

= B 

kl = 

'k 

0 

O' 

2 

3 2 



0 

k 

0 

.1 

2 1. 



.0 

0 

k. 


'(1 +k) 

1 

2 ' 





2 

(3 + k) 

2 





. 1 

2 

(1 +k)_ 
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Answers 


B—kl = 


3 P.Q = 


Q.R = 


P(RQ) = 


4 A.S = 


S.A = 


(1 -k) 1 

2 (3-fc) 

1 2 

8 4 0 
4 2 4 
11 5 7 

2 2 2 ' 

2 -3 -5 

10 8 12 ^ 

38 20 22 
12 6 4 

|_35 17 15 J 

6 6 4 

9 24 
-6 16 
4 2 

9 18 
-8 16 J 


—4 
2 
6 

-6 

4 


S 2 = 


2 
2 

(P0K = 


P(QR) = 


P(P0 = 


4.T = 


T.A = 


16 4 4 

12 10 14 
28 19 25. 

16 4 4 

12 10 14 
28 19 25 
-3 -1 -7 

53 25 25 

49 23 29 

6 -6 4 

—4 -9 -24 

2 6 -16. 
6 4 2 

6 -9 -18 
-4 8 -16 


"4 

0 

O' 

S.T = 

'4 

0 

O' 

0 

9 

0 


0 

-9 

0 

.0 

0 

16. 


.0 

0 

-16. 


5 A. V = 

"1 2 3' 

6 3-2 

V.A = 

' 1 

-2 

-2 4' 

3 6 


-1 

<N 

1 

_1 


. 3 

2 1. 


Examples XV 

1 |P| = -2 

2 iei = -6 











































































Introduction to Matrices 

3 |4|-5 

4 Singular 


5 C = A.B = 


6 8 

7 24 

8 |£>|-4 

Examples XVI 


8 13 13 

13 20 20 
7 14 14 


1 A' = 


2 B' = 


4 C = 


3 2 
2 3 
1 2 

-2 1 
3 -2 

6 4 

-4 -2 10 

1 2 -7 

2 -2 -2 


[adjoint Q\ = 


-4 

-2 


C is singular 


|£| = 4 


1 

2 


<?" = - 


5 = 5 
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10 

-7 

1 

'-4 

1 

2 ' 

6 

-2 

2 

-2 


. 10 

-7 

- 2 . 

r_ 

1 

3 

- 3 ' 


2 

-2 

-2 


\Q\ = -6 


-6 3 
2 -1 


V-. 


Answers 


6 P"‘ 



O' 

2 

- 4 . 


*“ = -8 


'-4 8 4 ' 
-6 4 2 
.-3 6 5 . 


A' = ' 

a 

h 

g’ 








h 

b 

/ 








Z 

f 

c. 







p Q = 


8 

5 

8" 

F = 

"l 

1 

3' 




4 

3 

4 


2 

0 

1 



- 

11 

8 

11. 


.1 

1 

2 . 


6' = 

"i 

2 

3' 


Q.P 


' 5 

3 

4 


i 

1 

2 




8 

6 

7 


j 

2 

3. 




14 

9 

11 

P Q -■ 


'5 

9 

14 

QP 

= 

'8 

4 

i r 



3 

6 

9 



5 

3 

8 



.4 

7 

11. 



.8 

4 

11. 


(P.0' = 

"8 

4 

ir 


5 

3 

8 


.8 

4 

11. 

(P0 = Q'P' 


and 

10 PQR = 

'26 

21 

47 


14 

11 

25 


.38 

30 

68 


(QPY = 


5 9 

3 6 

4 7 


14' 

9 

11 . 


(Qpy = P'0 




' 4 3 7 " 

6 5 11 
.10 8 18 . 
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Introduction to Matrices 


P' = 

*1 

1 

3" 


(QR)P' - 

'26 

14 

38' 


2 

0 

1 



21 

11 

30 


.1 

1 

2. 



.47 

25 

68. 

P'Q' = 


'4 

6 

10' 

RQ'P' = 

'26 

14 

38' 



3 

5 

8 


21 

11 

30 



.7 

11 

18. 


.47 

25 

68 j 


(P<2K)' = p6 14 38 
21 11 30 

.47 25 68. 
(PQR)' = (QR)'F = RQ'P' 

Examples XVII 


1 

'1 2 

r 

V 

= 

O' 


3 2 

i 

y 


2 


.2 3 

2. 

z_ 


.2. 



J 

p-‘ 

= 

1 

2 


x = 1 


y = 

-2 

1 

2 

'-1 

2 

O' 

fx" 

= 


-3 

1 

2 

y 



. 3 

0 

-4. 

|_z_ 



1 -1 


-r 

2 

-3. 
2 

3 

4 


-4 

5 

z = 3 
' 5 
-1 
9 


i.e. x = — 1 


1 -4 



"x" 


' 1" 


y 

= 

-2 


z, 


3. 


i.e. x = 2 


y = 2 
y = 3 


z = -3 


z = 4 
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zlnswws 


x" 

= 

V 



3 

_z_ 


.2. 


i.e. x = 1 y = 3 z = 2 


Examples XVIII 

1 


2 

3 

4 
6 


M is singular 
MN is singular 

in = i 


1 A.D = 


2 B.D = 


3 in = 

4 |G| = 

5 BT = 


1 

0 

3 

= 5 

1 

0 

2 

2 

1 

4 


3 

4 

5 

0 

1 

3 


5 

6 

7 

3 

0 

0 





0 

4 

0 

= 24 




0 

0 

2 





Ml 

= 

4 


in 

= 



= -6 


|l/| = 1 


7 Ml = |ri/| = \t\.\v\ = i 
Examples XIX 


'6 


3 


-4" 


D.A = 


6 - 

•2 - 

2 


3 


0 




3 

3 l 

.4 


-6 


0. 



- “ 

4 

4 1 

1 

8 

21' 



D.B = 

'1 

2 

3' 

2 

0 

28 




8 

0 

16 

,1 

8 

42. 




.7 

14 

42. 

-3 


in 

= 


6 

\EF\ = 

18 



-12 


|G 2 | 

= 

144 

|3G| 


-324 


'3 

2 

2' 




TB = 

'4 

0 

4' 

1 

2 

2 





1 

1 

2 

.3 

1 

0. 





.3 

-2 

0. 
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Introduction to Matrices 


6 TV = 


V 2 = 


7 A = LU 


-4 


|T| 

= 1 \BT\ = \TB\ 

= 

-4 

'2 

1 

O' 

VT = 

'1 

1 

O' 

1 

1 

0 


1 

2 

0 

.0 

0 

1. 


_0 

0 

1. 

'1 

0 

O' 

T 2 = 

"1 

2 

O' 

2 

1 

0 


0 

1 

0 

0 

0 

1. 


0 

0 

1. 


= 1 d e 

a (ad + 1) ( ae+f) 

.b (bd + c ) (be + cf+ 1)_ 

L and U are triangular matrices. 


8 

A = 

'1 

0 

O' 

'1 

2 

2' 



2 

1 

0 

0 

1 

4 



„3 

1 

1. 

.0 

0 

1. 

9 

L~ l - 


' 1 


0 

O' 





-2 


1 

0 





.-1 


-1 

1. 



10 |L| = 1 
|l/| = 1 


Examples XX 


1 



-4 

3 




2 

-2 3 

= 2 

Inverse = \ 

-7 -3' 

= 


4 -7 



.-4 -2. 
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Answers 


' 3 

-i 

O' 

or 

i 

3 

* 2 

-1 

O' 

0 

-1 

1 



0 

-3 

3 

--i 


- 1 . 



.-1 

5 

-3. 


5 ~ 


if-6 


6 -1 


-1 3 -3 

-6 3 2 

2 -1 1 
3 2 

2 -1 -2 
10 -7 —2J 
0 0 - 1 ' 
0-1 0 
-1 0 0 . 


or 


0 0 1 
0 1 0 
1 0 0 


Examples XXI 

1 \A-U\ = 0 

2 |B-A/| = 0 

3 A = 1 or -1 


M = I 1 0 
.0 3 
6 A = 4 or 3 


M = Q-'BQ 


A = 3 or 1 
A = 4 or 3 


4 Eigenvectors 

' 3 ' 

and 

2' 


-4 


. -3 . 


5 Q = 

3 2' 

Q' 1 = -1 

-3 -2' 


. 

1 

1 

1_ 


4 3. 


(2 = 

i r 

<r‘ = 

0 -1' 


. -1 0 . 


. 1 1 
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Introduction to Matrices 
’4 0 
.0 3. 

7 Since M = Q~'BQ , then M and B are similar matrices. 

Examples XXII 

1 Latent roots are 1, —1,2 


2 Eigenvectors are 

■ r 

■ r 

■ r 


-i 

0 

-l 


. 0 . 

,-i. 

_-i_ 


Q = 

l 

i 

r 

Q~ l = 

" i 

0 

r 


-l 

0 

-l 


i 

i 

0 


. 0 

-1 

-l. 


_-i 

-i 

-l. 


M = Q~ l AQ = [ 1 0 0 

0-10 
.0 0 2 . 

Latent roots are 2, 3, 4 


Eigenvectors are 


f 

-1 

0 


1 

0 

-1 


1 

-1 

-I 


Q = 

l 

l 

r 

Q~ l = 

' 1 0 

r 


-l 

0 

-l 


1 1 

0 


. 0 

-1 

-l. 


--1 -1 

-i. 


M = Q-'AQ = 


2 0 0 
0 3 0 
0 0 4 


Examples XXIII 


1 A(2 x 4) 

F(4 x 3) 


B( 2 x 2) 
G(3 x 3) 


C(2 x 4) D(4 x 2) E(4 x 4) 
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Answers 


2 A + C = [ 6 3 1 

L 2 10 9 

3 A.F (2x4).(4x3) = 

B.C (2 x 2H2 x4) = 
4 .D (2 x 4)(4 x 2 ) = 
E.F (4 x 4)(4 x 3) = 


F.G (4 x 3X3x3) = 


19 

5 

26 


(2x3) 

34 

10 

84 



21 

18 

21 

30 


. 7 

6 

7 

10. 


28 

5' 


(2x2) 

67 

63 . 




19 

5 

28' 


(4x3) 

38 

10 

62 



4 

-1 

10 



.27 

8 

46. 



10 

1 

18' 


(4x3) 

12 

-3 

30 



16 

-5 

30 



. 18 

-1 

38. 




C.D (2 x 4X4 x 2) 
A.E (2x4X4 x4) 


25 30 
112 43 
22 12 5 28 
16 26 5 35 


(2x2) 


(2x4) 


4 x' = x + 3 
/ = y + 4 

5 x' = x + 3 

/ = -y -4 


Translation 4 units upwards and 3 units hori¬ 
zontally. 

Translation as in question 4 followed by a 
reflection in the x-axis. 


Examples XXIV 

1 4x 2 + 6xy + 2y 2 

2 x 2 + 4xy+y 2 
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Introduction to Matrices 


[x y] 
[x y I 
[x y) 
[x y] 
[x y) 


1 1 
1 1 

2 3 


3 

2 

LI 

x' 

y. 

X 

y . 


= 25 


8 2x 2 + 5y 2 

9 x = — 3 or — ^ 

10 x 2 + y 2 +6x + 4v — 3 = 0 


11 

[x y 

1] 

'1 

0 

3 

X 




0 

1 

-2 

y 




.3 

-2 

9. 

.1. 

12 

[x y 

1] 

1 

0 

3' 

x' 




0 

1 

-2 

y' 




.3 

-2 

9. 

.1 j 


= 0 


= 0 


13 x 2 +y 2 +9x + 6y+3xz+2yz—3z 2 = 0 
Examples XXV 


1 A' = 

'2 4 

B' = 

-2 —3 


.3 7. 


2 2. 


MB)' = 

-13 

-29 

B'A' = 


. io 

22. 



2 C = 

" 2 3' 

(/IBCy = 


. -1 4. 



-13 -29 
10 22 
4 8 

53 117 
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Answers 


C(AB)' = 


' 4 

8' 

C’B'A 1 = 

' 4 

8' 

.53 

117 


.53 

117. 


B-'/r 1 = 


22 -10 
29 -13 


3 AB = 

-13 

10 

(AB)~ l =i 

'22 

-10" 


.-29 

22. 


.29 

-13. 


(C-‘) = 

1 

TT 

'4 

-3" 


(CT 1 


: TT 

'4 - 

-3' 




.1 

2. 





.1 

2. 


(a) P 1 = 

‘1 


2 O' 


Q' = 

'3 

4 

1* 




0 

1 1 



2 

2 

3 




.3 


4 3. 



.1 

2 

0. 



(PQ)’ = 

' 6 

14 

7' 

Q'P' ■■ 


' 6 

14 


7" 


11 

18 

11 



11 

18 

11 



1 

4 

2. 



. 1 

4 


2. 


(b) P-' =i 

'-1 3 —3 

-6 3 2 

Q~ l = -3 

'-6 3 2' 

2 -1 -2 


.2-1 1. 


. 10 -7 -2. 


(PQ) 


-1 _ _ 


ITT 


-8 -11 
0 
28 


26 
5 -10 
11 -46 


Q-'p-i 


= -JQ 


'-8 -11 26' 
0 5-10 

28 11 -46. 


7 (A — B)(A + B) = A 2 + AB — BA — B 2 and AB ^ BA 

8 I —A 2 . 
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Introduction to Matrices 


Examples XXVI 

1 Linearly independent 

2 Linearly independent 

3 No 

4 No. Linearly independent 

5 ad—be = 0 

6 p—2q + r = 0 A = 0 because they are linearly dependent 

7 No. Linearly independent 

8 Yes 

9 A = 0 when linear dependence is involved 

Examples XXVII 

1 2 

2 3 

3 2 

4 1 \D\ = 0 

5 2 |£| = 0 

6 4 


Index 
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This book is designed for the ordinary student 
rather than the university mathematician, and pre¬ 
sents the subject in a practical way. In the early 
stages mainly numerical methods are used to in¬ 
troduce the subject. First a study of simple (2 x 2) 
matrices leads to ideas of transformations or map¬ 
pings and then the concept of a vector is built up 
from numerical examples. The important ideas of 
Reflections, Stretches and Shears and their corres¬ 
ponding matrix forms are illustrated. The main con¬ 
cepts of (3 x 3) matrices are developed in a similar 
manner, including geometrical ideas in three di¬ 
mensions. Finally the student is introduced to some 
of the most important concepts of matrix algebra 
including Eigenvalues, Eigenvectors, Character¬ 
istic Equation, Linear Dependence and Rank. 
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